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1 INTRODUCTION 



ABSTRACT 

The Fermi-LAT collaboration has recently reported the detection of angular power above the 
photon noise level in the diffuse gamma-ray background between 1 and 50 GeV. Such signal 
can be used to constrain a possible contribution from Dark-Matter-induced photons. We es- 
timate the intensity and features of the angular power spectrum (APS) of this potential Dark 
Matter (DM) signal, for both decaying and annihilating DM candidates, by constructing tem- 
plate all-sky gamma-ray maps for the emission produced in the galactic halo and its substruc- 
tures, as well as in extragalactic (sub)halos. The DM distribution is given by state-of-the-art 
iV-body simulations of cosmic structure formation, namely Millennium-II for extragalactic 
(sub)halos, and Aquarius for the galactic halo and its subhalos. We use a hybrid method of 
extrapolation to account for (sub)structures that are below the resolution limit of the simula- 
tions, allowing us to estimate the total emission all the way down to the minimal self-bound 
halo mass. We describe in detail the features appearing in the APS of our template maps and 
we estimate the effect of various uncertainties such as the value of the minimal halo mass, 
the fraction of substructures hosted in a halo and the shape of the DM density profile. Our 
results indicate that the fluctuation APS of the DM-induced emission is of the same order as 
the Fermi-LAT APS, suggesting that one can constrain this hypothetical emission from the 
comparison with the measured anisotropy. We also quantify the uncertainties affecting our 
results, finding "theoretical error bands" spanning more than two orders of magnitude and 
dominated (for a given particle physics model) by our lack of knowledge of the abundance of 
low-mass (sub)halos. 



The isotropic gamma-ray background (IGRB) is the radiation that 
remains after the resolved sources (both extended and point-like) 
and the galactic foreground (produced by the interaction of cos- 
mic rays with the interstellar medium) are subtracted from the all- 
sky gamma-ray emission. A guaranteed component of the IGRB 
is the emission of unresolved known sources, whose contribution 
has been estimate d from population studies of their resolved coun- 
terparts: blazars ( Stecker et al. 1993; Stecker & Salamon 



Pavlidou & Ventersl2008l;llnoue & Totanil2009l;lAbdo et alfeOlOtt 
Abazaiian et alj201 ltlStecker & VentersfeOl lUsingal et al.ll2012h . 



star-forming galaxie s (Bhattacharva et al. 2009; Fields et al. 
Makiva et al] 1201 ll; lAckermann et all 1201 2b|; lLacki et all 
Chakrabortv & Fieldj|20l3T radio galaxies jstawarz et all 



2010; 



2012; 



2006; 



iMuecke & PohH 1 19981 ; iNarumoto & Totanil 12006^ 



Dcrmci 



1996; 



Massaro & Aiellol201ll; Inouel2011 ), pulsars and millisecond pul - 
sars dFaucher-Giguere & Loebl 20101; ISiegal-Gaskins et al. I l201lh . 
Gamma-Ray Bursts I Casanova et alj2008h and Type la Supernovae 
dLien & Fields! 1201 2E T Additional processes may also contribute 
to the IGRB such as cosmological structure formation shocks 



2007; 
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(e.g. lLoeb & Waxmanl l200d ; iGabici & Blasil 120031) . and interac- 
tions of cosmic rays (CRs) w ith the extragalactic background light 
(EBL) jKalashev et alj 120091) or with small solar system bodies 
dMoskalenko & Portej20(m 

Current estimates, however, suggest that the total unresolved 
emission from the classes listed above is not able to account for 
the whole IGRB intensity (e.g. lAielloll20T ]]), which strengthens 
the possibility that additional, unconfirmed sources are required to 
match the data. Gamma rays from Dark Matter (DM) annihilation 
or decay could explain the missing emission. 

DM is the dominant matter component of the Universe, re- 
sponsibl e for approximately one quarter of the energy density to- 
day (e.g. lJarosik et ali201 lh . We know little about its nature, apart 
from the fact that it has to be non baryonic. A well-studied class 
of DM candidates is that of Weakly Interacting Massive Particles 
(WIMPs), whose masses and interactions (set by the scale of weak 
interactions), offer promising non-gravitational signals for their de- 
tection in the near future. Within the context of annihilating DM, 
WIMPs are favoured by the fact that they naturally have a relic den- 
sity that matches the obser ved DM abundance (e.g. lKolb & Turner] 
1 1994 iBertone et ai1l2005l) . while for decaying DM, it has been 
shown that WIMPs can have a decay lifetime larger than the age 
of the Universe, and are therefore viab le DM candidates (see e.g. 
iBolz et al.ll200ll. lArvanitaki et alj|2009l) . WIMPs are also appeal- 
ing because their existence is predicted by fundamental theories 
beyond the Standard Model of Particle Physics, such as Super- 
symmetry (SUSY), Universal Extra-Dimensions or models with T- 
parity. In this paper we assume that DM is made of WIMPs, without 
making a specific assumption about the theoretical particle physics 
model from which WIMPs arise. 

This work is concerned with indirect detection of DM, i.e., the 
possibility of revealing the presence of DM from detection of its an- 
nihilation or decay products. In particular, we focus here on the case 
of gamma rays as by-products, studying the possible contribution 
to the IGRB coming from the DM annihilations (or decays) in the 
smooth DM halo of the Milky Way (MW) and its galactic subhalos, 
as well as from extragalactic (sub)halos. These contributions have 
already been es timated in the past using analytical and numerical 
techniques (e.g.lUllio et al.ll2002!;lTaylor & Silkll2003l;lAndd | 200l ; 
Ando & Komatsull200d; lAndo et al.ll2007bl: Isiegal-Gaskinsikood; 
Andol |2009l; iFornasa et alj 12009b IZavala et alj 120101: llbarra et al.1 



20091 : Isiegal-Gaskins et alj|201ll) and associated wi t h DM dAndo 
200a lAndo & Komatsul 120061: lAndo et al.1 |2007bl: ICuoco et al 
12001 l2008blTaoso et alj|2009; Siegal-Gaskins 2008; Fornasa et al 



20ld ; lHutsi et~al 2010; Cirelli et al. 201 1 ; Zavala et al. 201 1]). The 



recent Fermi-LAT measurement of the energy spectrum of the 
IGRB has been used to put constraints on the nature of the DM 
candidate by requiring that the DM-induced emission should not 
exceed the observed IGRB jAbdo et alj|2010al ; iHutsi et al.1 [20 id ; 
IZavala et alJl201ll ; ICalore et alj|2012h . The constraints derived are 
quite comp etitive: for instanc e, the most optimistic scenario con- 
sidered bv lAbdo et alj d2010al) puts a strong upper limit to the anni- 
hilation cross section, which is already of the order of the thermal 
relic value for a DM particle lighter than 200-300 GeV. 

The energy spectrum is not the only piece of information we 
can extract from the IGRB. Thanks to the good angular resolution 
of Fermi-LAT, it is also possi ble to measure its angular power spec- 
trum (APS) of anisotropies. lAckermann et alj d2012al) reported a 
detection of angular power in the multipole range between I = 155 
and 504 with a significance that goes from 7.2<x (in the energy bin 
between 2 and 5 GeV) to l.lcr (between 10 and 50 GeV), which 
represents the first detection of intrinsic anisotropies in the IGRB. 

There are different predictions for the normalization and 
shape of the APS produced by different populat i ons of unresolved 
sources, both asi 



2009; Siegal-Gaskins & Pavlidou 2009; Ando 2009; Zavala et al. 
bold ; llbarra etalj|20ld : ICuoco et al.1 1201 lh . The comparison of 



these predictions with the Fermi-LAT APS data can, in princi- 
ple, constrain the c ontribution of each source class to th e IGRB 
(ICuoco et alj|2012h . The analysis from lAckermann et alj d2012al) 
seems to suggest an interpretation in terms of a single population 
of unresolved, unclustered objects, due to the fact that the APS is 
roughly scale-independent over the energy range analysed. This 
recent measurement can then be used to complement other con- 
straints on a possible DM contribution to the IGRB. In the present 
paper, we take a first step in obtaining such constraints by revisiting 
and updating the prediction of the DM-induced emission (through 
decay and annihilation) and its associated APS, as well as estimat- 
ing the uncertainties involved. The comparison of these predictions 
with the Fermi-LAT APS data will be done in a follow-up study. 

In order to compute the DM-induced APS we combine the re- 
sults ofjw^J^body simulations of the galactic (Aquarius, hereafter 
AO. ISpringel et a l. 2008) and extragalac tic (Millennium-II, here- 
after MS-II. iBovlan-Kolchin et alj|2009l) DM structures, to con- 
struct all-sky maps of the gamma-ray emission coming from the 
annihilation and decay of DM in the Universe around us. Although 
we only focus here on the study of the anisotropy patterns in the 
gamma-ray emission, these maps represent per se a useful tool for 
future projects on indirect DM detection and we plan to make them 
available shortly after the publication of the follow-up paper dedi- 
cated to the comparison with the Fermi-LAT APS data. 

The extragal actic component i s expected to be almost 
isotropic (see, e.g.. IZavala et alj[20ld) , while the smooth galactic 
one is characterized by an intrinsic anisotropy, as a consequence 
of our position in the MW halo. The presence of galactic subha- 
los, however, reduces the expected gradient of the DM-induced 
gamma-ray flux as one moves away from the Galactic Center (GC). 
In fact, due to the large abundance of substructures and their more 
extended distribution, strong gamma-ray emission i s also expected 
quite far away from the GC (as it can be seen , e.g.. in | Springel et al 
2008blFornasa et alj2009l ; ICuesta et alJbOlU Sanchez-Cond e et al 



201 lh . 



Even though numerical simulations represent the most reli- 
able method to model the non-linear evolution of DM, they are 
limited by resolution. Since the minimum self-bound mass (M mm ) 
of DM halos is expected to be many orders of magnitude below 
the capabilities of current simulation!], this poses a challenge for 
an accurate predictio n and represen t s one of our largest sources 
of uncertainty (e.g Favlor & Silkl 120031; ISpringel et al.1 [20081; 
Siegal-Gasldnsl2008l;lAndcl2009l:lFornasa et al.l2009l:IZavala et al.1 

Sanchez-Conde et al.1 l201ll ; 
To address this problem, we 



20ld; iKam ionkowski et al. 



Pinzke et alj|201 lb iGao et aT] ; 



use a hybrid method that models the (sub)halo population below the 
mass resolution of the simulations by extrapolating the behaviour 
of the resolved structures in the MS-II and AQ simulations towards 
lower masses. Furthermore, we compute multiple sky maps with 
different values of M m j„ to determine with more precision what is 
the impact of this parameter on the the DM-induced emission. We 



1 The actual value of M m j n is related to the nature of the DM particle, 
with typical values coveri ng a quite large range, approximately between 



b produced by different populat i ons ot unresolved with typical values covering a quite large range, approximately be 

trophysical dAndo et alj|2007al ; lAndo & Pavlidod l(r 12 M and 1(T 3 M G (e.g. lProfumo et alj2006tlBringmannl2009h . 
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also consider possible effects due to different DM density profiles 
for the smooth halo of the MW. 

Such a detailed study of the uncertainties associated with the 
APS allows us to quantify, in addition to the normalization and 
shape of the APS, a "theoretical uncertainty band", that will prove 
to be useful in the comparison of our predictions with the Fermi- 
LAT APS data. 

The paper is organized as follows. In Sec. [2] we describe the 
mechanisms responsible for the gamma-ray emission from DM an- 
nihilation or decay. We then present how the data from the MS-II 
and AQ simulations are used to construct template maps of DM- 
induced gamma-ray emission from extragalactic DM (sub)halos 
(Sec.O and from the smooth galactic halo and its subhalos (Sec. 
[4^. In Sec.|5]we present the energy and angular power spectra, dis- 
cussing the different components and estimating their uncertainties. 
We discuss the implications of our results in Sec. [6] while Sec.|7]is 
devoted to a summary and our conclusions. 



2 DARK-MATTER-INDUCED GAMMA-RAY EMISSION 

In the case of DM annihilation, the gamma-ray intensity (defined 
as the number of photons collected by a detector per unit of area, 
time, solid angle and energy) produced in a direction *P is: 



dE 



&nm x 



p\A{z)^)e 



dN y (E y (\ + z)) 
dE 



(1) 



where E y is the observed photon energy, m x is the mass of the 
DM particle and (cr ann v) its annihilation cross section. The sum 
runs over all annihilation channels, each one characterized by a 
branching ratio, B', and photon spectrum (yield), dN' y /dE, com- 
puted at the energy of emission. The integration is over the line of 
sight (parameterized by A) to account for the redshift-dependent 
DM density field dA = c dzH(zY l . The exponential factor ac- 
counts for photon absorption from pair production due to inter- 
actions with the EBL along the line of sight, parametrized by an 
o ptical depth tebl(z, E y ), w hich we take from the model developed 
in Doming uez et alj {201 00. The first part of the integrand in Eq. 
[T]is usually referred to as the "particle physics factor" and only de- 
pends on the properties of DM as a particle, whereas the second 
part is called the "astrophysical factor" and depends on how DM is 
distributed in spac^E 

In the case of DM decay, Eq.Q]should be re-written as: 



d® 

-E {E ^ 



— - — f rf/lV B, 



dN' y (E y (l + z)) 
dE 



(2) 



p(Mz),V)e 



-TEBLfc.Ey) 



where the decay lifetime r is used instead of the annihilation cross 
section and the dependence on density is linear instead of quadratic. 
In the current section we describe the particle physics factor, 



introducing the mechanisms of gamma-ray production considered. 
Secs.[3]and|4]are devoted to the astrophysical factor. 

As mentioned in the Introduction, rather than considering a 
specific particle physics model, we focus on a general WIMP candi- 
date, which, for our purposes, is completely defined by m x , (cr aim v) 
or t, and its gamma-ray photon yield. The latter receives contribu- 
tions from three different mechanisms of emission: 

• Prompt emission: This radiation comes from the products 
of DM annihilation/decay directly, without any interaction with 
external particles. Within this first category, one can distinguish 
three different processes: i) gamma-ray lines from direct anni- 
hilation/decay into photons, ii) hadronization of quarks followed 
by neutral pion decay into photons and Hi) gamma rays from 
final state radiation and internal bremsstrahlung whenever there 
are charged final states. For DM annihilation, the branching ra- 
tios for monochromatic lines are usually subdominant and quite 
model-dependent (at least for SUSY models), while for DM de- 
cay, emission lines may be more prominent l Choi et al] |2010| : 
IVertongen & Wenigerl 1201 ll : iGomez- Vargas et alj I2012T) . In this 
work, we do not consider the emission from monochromatic lines, 
instead, we focus on mechanisms ii) and Hi), which are char- 
acterized by a continuum emission (e. g. |Fornengoetal.ll2004l : 



iBertone et al.l2005l ; lBergstrom et al.l2005klBringmann et alj2008l) . 
The continuum emission induced by hadronization shows some de- 
pendence on the DM mass and the particular annihilation/decay 
channel, but it is a mild one and the shape is more or less univer- 
sal. Finally, internal bremsstrahlung may also contribute inducing 
harder spectra and the possibility of bumps near the en ergy cut-off 
set by m x (see e.g. lBringmann et alj|2008tl2009l,l2012h . 



• Inverse Compton (IC) up-scattering: This secondary radiation 
originates when low energy background photons are up-scattered 
by the leptons produced by DM annihilation/decay. Since large y 
factors are required, usually one focuses on the case of electrons 
and positrons interacting with the Cosmic Microwave Background 
(CMB) photons and with starlight (either directly or re-scattered 
by dust). The amplitude of the IC emission and its energy spec- 
trum depends on the injection spectrum of e + / e~ and on the energy 
densi t y of the background rad i ation fields l Colafrancesco et all 
l2006llProfumo & Jeltemall2009l : Izavala et aljhoilh . For massive 
DM candidates, those IC photons can fall within the energy range 
detected by Fermi-LAT and, in some cases, represent a signifi- 
cant contribution to the DM-induced emission I Profumo & Jeltem3 
120091: iMeade et al.l20ld : iHutsi et alj201ol ; iPinzke et al.l201lh . See 



Appendix|A]for details on the computation of the IC emission. We 
note that for the case of extragalactic DM (sub)halos (Sec. [3} we 
only consider the CMB as a background source. This is mainly 
because the bulk of the emission comes from small (sub)halos 
(see Sec. 15. lb that are essentially em pty of stars and therefore 
lack any starlight background (see e.g. IProfumo & Jeltemall200^ : 
IZavalaetalJl201ll) . On the other hand, for the case of the MW 
smooth halo, a complete model for the MW radiation field is used 
(see Sec. 14. U. 



2 We have not checked the effect of other EBL attenuation models, since 
for the energies we consider in this work (from 0.5 GeV to 50 GeV), the 
contribution of the damping £~ t ebl factor is marginal. 

3 The particle physics and astrophysical factors are not completely inde- 
pendent from each other: the presence of M m j„, which is fixed by the parti- 
cle physics nature of the DM candidate, determines the minimum (sub)halo 
mass scale to be considered. Moreover, the dependence on redshift is both 
for the DM distribution and for the energy in the photon yield dN y /dE. 



• Hadronic emission: This radiation comes from the interac- 
tion of hadrons produced by DM annihilation/decay with the 
interstellar gas, and its contribution depends on the injection 
spectrum of hadrons and on the spatial distribution of ambi- 
ent gas. To implement this com ponent we follow the method 
described i n Delahav e et al. | d201lh (see also Icholis e7ai]|201ll : 
IVladimirov et al. 201 ll) and present the details of the calculation 
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Figure 1. Gamma-ray intensity from DM annihilation (solid lines) and de- 
cay (dashed lines) coming from the MW smooth halo (see Sec. I4.lt . For 
the "i-model" (black, red and yellow lines) the mass of the DM particle 
is 200 GeV for the case of annihilation and 2 TeV for decay. We assume 
(c<mnV) = 3 X 10~ 2s cm 3 s~' and r = 2 X 10 27 s, respectively. For the "r- 
model" (blue and green lines) the parameters are the same except for the 
mass which is 2 TeV for both annihilating and decaying DM. Black and 
blue lines indicate prompt-emission, red and green IC emission, and yellow 
hadronic emission. The latter is not shown for the r-model. 



in Appendix|B] We only consider this additional component for the 
case of the MW smooth DM halo. 



As benchmarks, in the remainder of this paper, we consider two 
commonly-used annihilation/decay channels, with which we illus- 
trate the role of the different mechanisms: a "i-model" for annihi- 
lation/decay entirely into bb quarks (B b =1) and a "r-model" for 
annihilation/decay into t + t~~ (B r = 1). The photon and e + Je~ yields 
are computed using the tables presented in lCirelli et al.l d201ll) . For 
both cases, we fix the annihilation cross section and decay life time 
to 3 x 10~ 26 em 3 s -1 and 2 x 10 27 s, respectively. The DM mass is 
selected to be 200 GeV for the b channel in the case of annihilating 
DM and 2 TeV otherwise. These values are chosen to be slightly 
below the most recent exclusion limits set by the Fermi-LAT d ata 
dAckermann et al.ll201 lkfougger et al.ll20ld : lHuang et alj|201 lh . 

In Fig.Q~]we compare the gamma-ray production mechanisms 
listed above. The lines indicate the energy spectrum of the emis- 
sion from annihilation (solid) or decay (dashed) of DM in the 
MW smooth halo (see Sec. 14. It . For the fo-model, prompt emission 
(black lines) always dominates over IC (red lines) and hadronic 
emission (yellow lines), both for annihilation and decay. On the 
other hand, for the r-model, IC (blue lines) overcomes the prompt- 
emission (green lines) at low energies. For the r-model, hadronic 
emission is negligible and is not plotted. 



3 THE GAMMA-RAY EMISSION FROM 
EXTRAGALACTIC (SUB)HALOS 

3.1 Resolved (sub)halos in the Millennium-II simulation 
(EG-MSII) 

The MS-II follows the formation and evolution of DM structures in 
a comoving cube of L = 100 Mpc//i on a side and a total of (2160) 3 
simulation particles (Bovlan-Kol chin et al.l|2009|) . The simulation 
is done within the context of a WMAP1 cosmology with the fol- 
lowing parameters: £l m = 0.25, fl A = 0.75, h = 0.73, cr 8 = 0.9 and 
n s = 1 ; where Q.„, and Q A are the contribution from matter and cos- 
mological constant to the mass/energy density of the Universe, re- 
spectively, h is the dimensionless Hubble constant parameter at red- 
shift zero, n s is the spectral index of the primordial power spectrum, 
and <x 8 is the tins amplitude of linear mass fluctuations in 8 Mpc//i 
spheres at redshift zero. Its mass resolution is 6.89 x l0 6 M B /h and 
there are 68 snapshots recording the particle distribution at different 
redshifts between z = 127 and z = 0. 

Instead of working directly with the particles in the simula- 
tions, we use the MS-II (sub)halo catalogs, which are const ructed 
using a friend-of-fri ends (FOF) algorith m toavis et"al]|l985l) and 
the SUBFIND code JSpringel et al.l200ll) that identifies self-bound 
substructures within FOF halos. Dealing with the (sub)halo cata- 
logs, instead of the particle data, has two advantages: it is much less 
expensive computationally and, more importantly, it avoids resolu- 
tion effects near the centre of DM (sub)halos, where the simulation 
particles severely underestimate the DM density (note that this is 
precisely the region with the highest gamma-ray production rate). 
On the other hand, we are neglecting the contribution from the DM 
mass that does not belong to (sub)halos. The emission rate from 
unclustered regions, however, is likely to be negligibl e especially 
for DM annihilations (see, e.g. Angulo & White 2009 who analyt- 
ically estimated that between 80-95% of the mass is in collapsed 
objects. In the case of decaying DM this suggests that, by neglect- 
ing unbound particles, we underestimate the luminosity by, at most, 
20%). 

The MS-II (sub)halo catalog contains the global properties 
needed for each object: its virial mass Mooo (defined as the mass 
up to rooo, w h ere the enclosed density is 200 times the critical den- 
sity), its maximum circular velocity V max and the radius r max where 
this velocity is attained. The latter two quantities completely deter- 
mine the annihilation/decay luminosity for each halo if we assume 
that they have a spherically symmetric d ensity distribution giv en by 
a Navarro-Frenk- White (NFW) profile jNavarro et alll 19971) . The 
number of gamma rays (per unit of time and energy) coming from 
a (sub)halo with a boundary at r 2 oo is then given by L = /ppL', 
where: 



/pp - 



L' 



(0-annV) 



2m 2 



Yb^ 



V 4 

1.23- max 



G^r ma x 



1 - 



1 



(1 + C 2 00) 3 



for the case of annihilation, and 
m v T t—t dE 



L — -Ldecay — 2.14- 



for the case of decay. 



ln(l + c 20 o) ■ 



''200 



1 + C200 



(3) 



(4) 



(5) 



(6) 
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The concentration C200 is a ls° determined from V ma x and /- max 
inverting the following relation (e.g. lSpringel et al . 2008): 

14.42 6 f-^U 2 = 292 C ™ . (7) 

\H(z) r max ) 3 ln(l + c 20 o) - c'20o/(l + c 2 oo) 

The choice of the NFW density profile is motivated by its uni- 
versality and by the fact that it gives a good fit to simulated DM 
(sub)halos over a large mass range. However, assuming other DM 
density profiles could have an impact on the total gamma-ray emis- 
sion, as well as on the shape and normalization of the APS. A dis- 
cussion about this possible source of uncertainty is left for Sec. [6] 

We define as EG-MSII the signal coming from (sub)halos in 
the MS-II catalogs with at least 100 particles; below this num- 
ber, the mass and abundance of DM objects in the MS-II can be- 
come unreliable. This sets an "effective" mass resolution of M lcs = 
6.89 x 10 8 M G //? for the extragalactic contribution. DM structures 
with less than a few thousand particles can be affected by numer- 
ical effects_(gj^ryjtational softening and two-body relaxation, see 
e.g. iDiemand et aT]|2004) that could influence the values of V raax 
and r raax and, as a conse quence, L ann or Ld eC a v We implement the 
prescription described in IZavala et al. to correct for these 

effects. 

In order to simulate the past light cone we need to probe a vol- 
ume which is much larger than the MS-II box. To do t his, we fol- 
low closely the procedure given in IZavala etaljfeoich which can 
be summarized as follows. The region around the observer is di- 
vided into concentric shells, each of them centered in redshift space 
on the discrete values Zt corresponding to each simulation output. 
The volume defined by each shell has a fixed size in redshift space 
and a corresponding comoving thickness which is filled with iden- 
tical, no n-overlapping copie s of the MS-II box at the redshift zt (see 
Fig. 9 of IZavala et alJl201Ql) . In order to compute the DM-induced 
gamma-ray emission from a given direction x ¥, we follow the line of 
sight defined by "P that crosses the MS-II replicas, and sum up the 
emission produced in all the (sub)halos encountered. The projec- 
tion into a 2-dimens ional map is done with the HEALPix packag^] 
dGorski et al.ll2005l) . assuming N_side=512, corresponding to an 
angular area of approximately 4 x 10 -6 sr for each pixel. If a given 
halo subtends an area larger than this value, then it is considered as 
an extended source. In this case, each of the pixels covered by the 
particular halo is filled with a fraction of the total halo luminosity, 
assuming the corresponding projected surface density profile. 

To avoid the repetition of the same structures along the line 
of sight (which would introd uce spurious periodicity along this di- 
rection), IZavala et alj d2010h used an independent random rotation 
and translation of the pattern of boxes that tessellates each shell. 
This method, however, still leaves a spurious angular correlation 
at a scale AO corresponding to the comoving size of the simula- 
tion box, which mainly manifests itself as a peak in the APS cen- 
tered on £* = 2n/A6. This angular scale decreases as we go deeper 
in redshift, since each copy of the MS-II cube covers smaller and 
smaller angles. This implies that the periodicity-induced peak will 
be located at a different multipole for each shell. Once the contri- 
butions from all shells are added up, this effect is largely averaged 
ou t, and the total APS is free from any evident features (see Fig. 12 
of IZavala et al.ll201oh . Nevertheless, the spurious angular period- 
icity, in addition to the fundamental angular correlation associated 
with Ad, introduces smaller scale harmonics that affect multipoles 
larger than £* . Although these additional peaks are much smaller 




Multipole 

Figure 2. Fluctuation APS of anisotropics of the gamma-ray intensity 
produced by DM annihilating in extragalactic (sub)halos resolved in the 
MS-II, and located in a shell corresponding to z < 0.01 (black line), 
0.19 < z < 0.22 (red and yellow lines) and 1.57 < z < 1 .70 (green and blue 
li nes). Yellow and b lue lines refer to the map-making algorithm presented 
in lZavala et alj feplOh . while the black, red and green lines correspond to 
our improved algorithm; see text for details. The grey band around the black 
line indicates the l<x standard deviation among 10 different realizations of 
the first shell (z < 0.01). 

than the fundamental one, we decided to reduce this spurious effect 
by randomly rotating and translating every single replica within the 
past light cone instead of doing so only for every concentric shell. 

The improvement of the new method becomes evident in Fig. 
[2]where we show the comparison between the fluctuation APs|f| (for 
individual shells) computed with ou r map-making code (red and 
green lines) and the original one by IZavala et alj fl2010l) (yellow 
and blue lines). The yellow and red lines refer to the shell with 
z = 0.21, while the blue and green lines are for z = 1.63. The small- 
scale spurious harmonics essentially disappear in the new method 
and the fundamental mode, although still present, is greatly reduced 
relative to the previous method. 

The map-making code produces realizations of the distribu- 
tion of DM halos around the observer through random rotations 
and translations of the MS-II boxes that fill the volume of the past- 
light cone. In order to quantify the effect of this random component 
in the simulated signal, we generate 10 different realizations of the 
first shell (corresponding to z < 0.01) and compute the fluctua- 
tion APS for each of them. We only consider the effect of having 
different random rotations for the first shell since it is expected to 
be more important for nearby resolved structures, while shells at 
larger redshifts are less affected. In Fig. [2] we plot the average APS 
over these 10 realizations (black line) as well as the l<x fluctuation 
(grey band). We can see that the effect induced on the APS is rela- 
tively small (at least compared to the other sources of uncertainties 
introduced later) and we neglect it from now on. 

All halos up to z = 2 are considered when computing the ex- 
tragalactic signal. By this redshift, the cumulative emission has al- 
ready reached > 80% of the total signal (in the case of prompt 



4 http://healpix.jpl.nasa.gov 5 The APS will be formally defined in Sec. 15. 21 
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emission ~ 90% of the signal actu ally comes from z < 1; see 
Fig. 9 of IProfumo & Jeltemal 120091 and Fig. 1 1 of IZavala et al.l 
|2010|) . The first shell of the extragalactic map starts at a distance 
of ^min = 583 kpc, corresponding to approximately twice the virial 
radius of the galactic halo. The volume within this distance is filled 
with the data from the AQ simulation (see Sec.|4](. 

In the upper panels of Fig. [3] we show the gamma-ray inten- 
sity of the EG-MSII component at 4 GeV for the first snapshot (z < 
0.01), in the case of annihilating DM (left panels, m x = 200 GeV, 
((TanV) = 3 x 10~ 26 cmV' and B h = 1) and decaying DM (right 
panels, m x = 2 TeV, t = 2 x 10 27 s and Sj, = 1). The characteristic 
filaments of the cosmic web and individual DM halos are clearly 
visible, as well as (at least for the case of decaying DM) some sub- 
halos hosted in large DM clumps. In the second row of this figure, 
we show the intensity up to z = 2: the map is much more isotropic, 
even if some of the prominent, closest structures can still be seen. 



3.2 Unresolved main halos (EG-UNRESMain) 

We describe now how we model the contribution of unresolved 
main halos, (i.e. those with masses below M lcs ), a contribution that 
we call EG-UNRESMain. Since this is a regime which goes below 
the MS-II mass resolution, we are forced to resort to some assump- 
tions concerning both the distribution and the individu al properties 
of DM halos. Our approach is similar to the one of IZavala et al.l 
(2010): we use main halos in the MS-II to perform an analytic fit 
to the single-halo luminosity (i.e. L(M) defined in Eqs.[4]and[6} as 
well as to the following function: 



F(M) ■■ 



ZUM) 
MA log M 



« In 10L(M) 



An(M) 
AM ' 



(8) 



which is the total luminosity of main halos with a mass in the log- 
arithmic mass range logM + AlogM/2, divided by its mean value 
M and the width of the logarithmic mass bin. The second equal- 
ity shows how F(M) depends on the halo mass function An/ AM 
in the bin considered. By extrapolating the fit obtained for F(M) 
above M rcs , it is possible to estimate the gamma-ray intensity due 
to DM main halo s below M KS down to different values of M min . In 
IZavalaetalJfeOld) (see their Figs. 5 and 6), the authors compared 
the predictions of such an extrapolation, in the case of annihilation, 
with the result of an analytical model based on the formalism of 
IShethetalJl200lh for the halo mass function an d lEkeetal .1 (2001) 
for the concentration-mass relation. The total flux in unresolved 
main halos with a mass between M m j n = 10~ 6 M G /7i and M res agrees 
within a factor of 5 between the two approaches. 

This missing flux is then added to the emission of resolved ha- 
los with mass between 1.39 x 10 s and 6.89 x l0 9 M e /h (halos with 
particle number between 20 and 100). The decision to boost up only 
these halos is equivalent to assuming that halos smaller than M lcs 
share the same spatial clustering as those within that mass range. 
This assumption is motivated by the fact that the two-point corre- 
lation function of halos app roaches an asymptotic valu e already at 
these masses (see Fig. 7 of lBovlan-Kolchin et al]2009T) . 

Nevertheless, the actual clustering of low-mass main halos 
is unknown and, even if they trace the distribution of more mas- 
sive objects, treating their contribution simply as a boost factor for 
the halos with lowest masses in the MS-II may overestimate their 
true clustering. Assuming, instead, that they are distributed more 
isotropically would reduce their contribution to the total APS (es- 
pecially at low multipoles), although it is difficult to estimate pre- 
cisely by how much. In what follows we assume that the uncer- 



tainty in the clustering of unresolved main halos is small and can 
be ignored. 

The same procedure described above is also applied for the 
case of DM decay. 



3.3 Unresolved subhalos 

In this section we describe how we account for the emission from 
unresolved subhalos, i.e.: i) subhalos with masses below M lcs that 
are hosted by main halos in the MS-II catalogs, and ii) subha- 
los of unresolved main halos. We do not consider sub-subhalos 
since their contribution is likely negligible in comparison (see, 
e.g., iMartinez et al. 2009). Note also that, at least at low red- 
shifts, most of these subhalos have been removed by tidal stripping 
dSpringel et alj2008l) . 

For the extragalactic emission, the impact of unresolved sub- 
halos on the intensity and angular anisotropy spectra is essentially 
to boost the luminosity of the host by a certain amount. Thus, in 
principle, any method that provides boosts within the range of what 
has been found previously in the literature is a reasonable one. The 
method we use here has the advantage of having a single parame- 
ter (k in Eq.[9]l that controls the abundance of substructure, which 
can be easily adjusted to obtain the subhalo boosts that have been 
repo rted in the past. 

iKamionkowski & Ko ushiappas d2008l) and 

iKamionkowski etlri] ^20 id) propose a method to compute 
the subhalo boost factor for the annihilation rate of a MW-like DM 
halo, providing an expression for the total boost factor B alm (M MW ), 
as well as the differential profile B ann (Mviw, r) (expressing the 
boost factor at a distance r from the center of the halo). This 
prescription was cal ibrated with the Via Lactea II simulation 
jDiemand et alj |2008). The distribution of particles in this simula- 
tion is used to derive the probability P(p, r) of having a value of the 
DM density between p and p + dp at a distance r from the center 
of the main halo. Two different components contribute to Pip, r): 
the first one is Gaussian and corresponds to the smooth DM halo, 
while for higher values of p, the probability is characterized by 
a power-law tail due to th e presence of subhalos (see Fig. 1 of 
IKamionkowski et alj2O10h . The fraction of the halo volume that is 
filled with substructures is well fitted b}[|: 



p sm (r) 



p sm (r = 100 kpc) 



(9) 



with k = 7 x 10~ 3 . P{p, r) is then used to derive an expression for 
the boost factor B mm (r) in the case of annihilating DM: B ann (M, r) = 
fT™* dpP(p, r)p 2 /p 2 m (r), where p max is a maximum density, which 
is of the order of the density of the earliest collapsing subhalos (see 
below ) and p ml is the density of t he smooth component. 

Sanchez-Conde et "aH d201lh extended the previous method to 
halos of all sizes, adopting a slight modification to the definition of 
Mr): 



l-/,(r) = * 



p sm (r = 3.56 x r s ) 



(10) 



Since current simulations are many orders of magnitude away from re- 
solving the whole subhalo population down to M m [ n , f s (r) is known with 
limited precision and represents one of the implicit uncertainties of our pre- 
dictions. 
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Figure 3. All-sky maps of the gamma-ray intensity (in units of cm s sr GeV ) at 4 GeV from DM annihilation (left panels) and DM decay (right panels). 
The figure shows the emission of all DM (sub)halos down to the resolution limit of the MS-II (EG-MSII component). In the upper row only nearby structures 
(z < 0.01) are considered, while in the second row the emission up to z = 2 is considered. In the last row we plot the emission from all extragalactic (sub)halos 
(resolved and unresolved) down to = 10~ 6 M o //i with the LOW subhalo boost (see text for details). In all cases, annihilation or decay into b quarks is 

assumed: for annihilating DM, m x = 200 GeV with a cross section of 3 X 10~ 26 cm 3 s~' , while for decaying DM, m x = 2 TeV with a lifetime of 2 X 10 27 s. The 
photon yield receives contributions from prompt emission and IC off the CMB photons (see Sec. [2}- In each map we subtract the all-sky average intensity of 
that component, after moving to a logarithmic scale. Note the different scales in the first row. 



where r s is the scale radius of the host halo given in We note 
that this implies that halos of all masses have the same radial de- 
pendence of f s , only rescaling it to the particular size of the halo. 
This is partially supported by the mass-inde pendent radial distri - 
bution of subhal os found in simulations (e .g. lAngulo et al]|2008h . 
Using Eq.[T0] S anchez-Conde et alj l l201ll) found that 5 aim < 2 for 
the MW dwarf spheroidals, while B„ ~ 30 - 60 for galaxy clus- 
ters (integrating up to the tidal and virial radius, respectively). In 
both cases, the morphology of the total gamma-ray emission com- 
ing from the halo is modified since the subhalo contribution makes 
the brightness profile flatter and more extended. 

For the case of annihilating DM, we account for the contri- 
bution of unresolved subha los by implementing the procedure of 
Sanchez-Cond e~etal1 ^201 ll) in two different ways: 



7 The value of 3.56 is chosen so that, for the MW halo in Via Lactea II, 
Eqs.[K)]and[9]are identical. 



• for the subhalos of unresolved main halos we integrate 
F m „(M)B mm (M) to compute the total luminosity from M mm to M m!i . 
The result of this integral is then used to boost up the emission 
of main halos in the MS-II with a mass between 1.39 x 10 8 and 
6.89 x 10 9 M e /h. 

• for subhalos belonging to main halos that are resolved in the 
simulation we boost up the luminosity of each halo by the mass- 
dependent boost S ann (M) (i.e. the integral of B ann (M, r) up to the 
virial radius). If the halo is extended, in addition to a total lumi- 
nosity boost, we assume a surface brightness profile as given by 
BsamiM, r). We need to apply a correction to this procedure since 
these equations account for subhalos from a minimum mass M mm 
up to the mass of the main halo M, whereas subhalos with masses 
above M lcs are resolved and already accounted for in the simula- 
tion (they belong to the EG-MSII component). To correct for this 
double-counting, we simply compute (and subtract) the emission 
due to subhalos down to a minimal mass equal to M m j„ = M ms . 
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We note that changing M mi n corresponds to changing p max . 
Fromfkamionko wski et alj fcoid) . the maximum density in a halo 
is the density that its smallest subhalo had at the moment this sub- 
halo formed: 



/ max 



200 c| 00 (M mi 
12 f(c 2 oo(M mm ,z F )) 



Pcritfe), 



(11) 



where /(x) = ln(l + x) - x/(l + x). The epoch of collapse zf 
as a function of halo mass can be computed using the spherical 
collapse model of DM ha lo formation and evolution (see, e.g., 
ISanchez-Conde et al.ll2007l and references therein), which shows 
that for low masses, up to ~ 1 M e /h, all halos collapse approxi- 
mately at the same redshift, zf = 40. The initial concentrations are 
set by the formation epoch, which means that halos that collapse at 
roughly the same zf will have similar c 2 oo(zf)- Thus, according to 
Eq.QT] all low-mass subhalos will be characterized roughly by the 
same p max ~ 2.51 x 10 9 M o /kpc 3 (for M min < lM e /h), after fixing 
C2qo(zf) to a constant value of 3.5 as su ggested by simulations (e.g. 
iDiemand et alj2006llZhao et al.l2 009 n. We compute p max for a set 
of reference values of M mm (see also Sec. 16.41 noting that zf > 5 
for M mm < l 9 M n /h, which im plied that we can safely assume 
C2oo(z F ) = 3.5 dZhao et al.f 2009). Note also that, by using the case 
of Mmin = M res = 6.89 x 10 9 M e /h we can correct for the aforemen- 
tioned problem of double-counting the subhalos with masses above 

the MS-II mas s resolution. 

Recently, IPinzke et all d201lh and iGao et al.l d201lh also es- 
timated the substructure boost for DM halos of mass ranging 
from those of dwarf spheroidals to those of galaxy clusters. They 
point to substantially large r boost factors than those found by 
Sanch ez^Conde et"afl d201lf) for the same mass range. This is 
mainly a consequence of the different methodologies. In the former 
cases, the subhalo mass function and the concentration-mass rela- 
tion are power-laws calibrated at the resolved masses and extrapo- 
lated to lower unresolved m asses. On the contrary, in the method by 
Kamionkowski et al. ( 2010) (with the modification implemented in 
Sanchez-Conde et"al]|201ll) . the dependence on M min is flatter to- 



wards lower masses due to the limit on the natal concentrations. 

Nevertheless, using the procedure described in the previous 
pa ragraphs, we can ob tain similar subhalo b oosts to those given 
bv lPinzkeetalJ ( l20TTh and IGao et al] J20l"ih if we substantially 
increase the parameter that controls the abundance of substruc- 
ture in Eq. [TOl to k = 0. 15. Both cases (k = 7 x 10~ 3 as in 



ISanchez-Conde et alj|2011, and k = 0.1 5 to reproduce the results 
of lPinzke et alj201 ll and lGao et alj201 l|) are considered in this pa- 
per as representative of scenarios with a small and a large subhalo 
boost and are referred in the following as the LOW and HIGH sce- 
narios, respectively. These two cases represent the extreme values 
reported in the literature for the contribution of unresolved sub- 
halos. By obtaining predictions for the total DM-induced emission 
for these extrema, we aim at estimating how large is the uncertainty 
associated with the unresolved subhalo population. Parameterizing 
such uncertainty in this way represents a "hybrid" approach, since 
it does not rely comple tely either on a d irect extrapolation of the re- 
sults of simulations (Zavala et al. 20 lot) or on analytical est imates 
such as the stable clustering hypothesis dAfshordi et alj20ld> . 

Up to now, the discussion of how to model unresolved subha- 
los refers only to the case of annihilating DM. For decaying DM, 



Here, a matter power spectrum parametrized as in iBardeen et all jl986h 
was used to compute zf , with the most recent values of the cosmological 
parameters and with no exponential cut-off at the minimal mass of DM 
halos. 



there is no need to model this contribution since these subhalos are 
too small to be detected by the subhalo finder and their mass is 
already accounted for in the mass of the host halo. Since for decay- 
ing DM the total luminosity of a halo is proportional to its mass, the 
unresolved subhalos contribute to what we call the "smooth com- 
ponent'^]. This is strictly valid only if we consider the total halo 
luminosity. If the intensity profile is needed, we should consider 
that the true spatial distribution of unresolved subhalos is expected 
to be different from that of the smooth component. In the case of 
the extragalactic emission we neglect this effect since only the halos 
that are close by appear extended in the maps, while the vast ma- 
jority appear as point sources. For the case of the galactic emission 
we comment on this issue on Sec. 14.21 



4 THE GAMMA-RAY EMISSION FROM THE 
MILKY- WAY HALO 

4.1 The smooth Milky-Way halo 

Our model for the emission from the smooth DM halo of our own 
alaxy is partially based on the results of the Aquarius project 
Springel et ai]|2008l : lNavarro et alj|2008h . With the goal of study- 
ing the evolution and structure of MW-size halos, the Aquarius 
project selected a group of MS-II halos with properties similar to 
the MW halo and resimulated them at increasing levels of reso- 
lution. The different AQ halos are characterized by virial masses 
between 0.95 and 2.2 x l0 12 M o /h and have a variety of mass accre- 
tion histories jBovlan-Kolchin et al.ll201oh . In this sense, they are 
not expected to be a perfect match to the dynamical properties of 
our own MW halo, but rather to be a representative sample of MW- 
size halos within the context of the CDM paradigm. We consider 
here the halo dubbed Aq-A- 1 , containing more than one billion par- 
ticles within T2oo and having a mass resolution of 1250 M Q /h. A 
careful analysis of the density profile of the s moo th component of 
the Aq-A-1 halo performed bv lNavarro et al.l d2008h shows that the 
simulation data is best fitted by an Einasto profile (preferred over 
an NFW profile): 



In 



pW 

P-2 



— I -1 

r-2 



(12) 



with r_ 2 ~ 15.14 kpc,p_ 2 = 3.98 x 10 6 M o /kpc 3 and a ~ 0.170. 

Stellar dynamics and microlensing observations can be used 
to constrain the absolute value of the DM density at the position 
of the Earth, p ]lx . Different result s point towards a range of values 
between 0.2 and 0.5 GeV/cm 3 dPrada et alj|2004l: ICatena & uliiol 
l20ld: IPato et alj|20ld : ISalucci et alfeOloUlocco et alj|201lh . NoT 
ing that a different value for the local DM density would shift up or 
down our predictions for the intensity of the emission from DM an- 
nihilation/decay in the MW smooth halo proportionally to pf oc and 
Pioc, respectively, we decide to renormalize the value of p_ 2 of Aq- 
A-1 in order to reproduce a refer ence value of p \ oc = 0.3 GeV/cm 3 
(a similar approach was used in lPieri et alj|201 ll) . 

To build our template map for the smooth MW halo, we as- 
sume that the observer is located at the solar circle at a distance of 
8.5 kpc from the GC and we integrate the DM-induced emission 



9 For the case of DM annihilation note that, although the mass of unre- 
solved subhalos is also accounted for as part of the "smooth component", 
this does not imply that their contribution to the gamma-ray intensity is al- 
ready considered since the annihilation rate is not proportional to the DM 
density, but to the density squared. 
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Figure 4. All-sky map of the galactic gamma-ray intensity (in units of cm 8 sr -1 GeV ) at 4 GeV from DM annihilation (left panels) and decay (right 
panels). In the first row, we show the emission from the smooth MW halo, while the contribution of resolved subhalos in the Aquarius Aq-A-1 halo (GAL-AQ 
component) is shown in the second row. The maps on the last row indicate the total galactic emission accounting for the MW smooth halo and its (resolved and 
unresolved) subhalos down to M ra j n = 10~ 6 M Q //i (for the LOW subhalo boost). As in Fig. [3] m x = 200 GeV, the cross section is 3 X 10~ 26 cm 3 s~' and Bj = 1 
for the left panels, while m x = 2 TeV with a lifetime of 2 X 10 27 s and B\, = 1 for the right ones. The intensity includes contributions from prompt emission and 
IC with the CMB photons (see Sec. [2). For the emission of the MW smooth halo we also consider IC with the complete ISRF, as well as hadronic emission. 
The non-prompt emission alone is shown in the smaller panels overlapping with the maps of the first row. In each map we subtract the all-sky average intensity 
of that component, after moving to a logarithmic scale. Note the different scale in the different panels. 



along the line of sight up to a distance of 583 kpc (~ 2.5 r 2 oo of Aq- 
A-1). This distance marks the transition between our galactic and 
extragalactic regimes and it is selected because the Aq-A-1 halo is 
still simulated with high resolution up to this radius, and it there- 
fore provides a better representation of the outermost region of the 
MW halo than the MS-II. For the smooth component, in addition 
to the prompt emission and secondary emission from IC scatter- 
ing with the CMB photons, we also consider the emission due to 
IC scatterin g with the complete Inter stellar Radiation Field (ISRF) 
provided in Moskalen ko et al.l d2006h as well as hadronic emission 
from interactions with the interstellar gas (see Appendices [A] and 
[B]for details). The first row in Fig. [4] shows the gamma-ray emis- 
sion from DM annihilation (left panel) and decay (right panel) in 
the smooth MW halo. The secondary emission correlated with the 
MW ISRF and the interstellar gas can be seen along the galactic 
plane and is plotted independently in the small panels overlapping 
with the maps of the first row. 



4.2 The Milky Way subhalos (GAL-AQ and GAL-UNRES) 

This section focuses on the contribution of galactic subhalos, deal- 
ing with i) subhalos that are resolved in the Aq-A-1 halo, (which we 
refer to as the GAL-AQ component) and ii) subhalos with masses 
below the mass resolution of AQ (which we call the GAL-UNRES 
component). As we did in Sec. 13. II we use the subhalo catalog to 
compute the luminosity of each object from its V max and r max val- 



ued I. Only subhalos with more than 100 particles are considered, 
resulting in an "effective" AQ mass resolution of 1.71 x 1O 5 M . 
The gamma-ray intensity in a given direction Y is then obtained by 
summing up the contribution from all subhalos encountered along 
the line of sight, up to a distance of 583 kpc. The GAL-AQ compo- 



10 As in the case of extragalactic (sub)halos, we correct the values of V max 
and r max for numerical effects (see Sec. l3.lt . 
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nent is shown in the second row of Fig.|4]in the case of annihilation 
(left) and decay (right). 

For an annihilating DM candidate, the contribution of unre- 
solved galactic subhalos is accounted for using the same procedure 
as for unresolved extragalactic subhalos described in Sec. 13.31 in- 
troducing the LOW and HIGH cases as representatives of scenar- 
ios with a small and a large subh alo annihilation boost . The L OW 
boost is taken again directly from lKamionkowski et alj d2010h and 
ISanchez-Conde et al.l d201ll) (which assumes k = 7 x 1(T 3 ), while 
the HI GH boost is tuned to reproduce the results of Spring el et al.l 
( 2008) who estimated a total subhalo boost of 232 (integrating up to 
r2oo for the Aq-A-1 halo and for M ra j n = 10~ 6 M o /A); we reproduce 
this result using k = 0. 

In the case of decaying DM, we note that the mass contained 
in resolved subhalos is 2.7 x 10 U M G (~ 15% of M 2 oo for the Aq-A- 
1 halo, if we conside r subhalos down to 1.71 x 10 5 M o , see Eq. 5 of 
ISpringel etall d2008h . This goes up to 3.9 x 10 1 1 M e if we extrapo- 
late the subhalo mass function down to M min = 10~ 6 M Q //i, which 
implies that unresolved subhalos contribute to the halo mass (and 
hence to the total decay luminosity) slightly less than resolved ones 
(see end of Sec. 13.31 ). Thus, an upper limit to the gamma-ray inten- 
sity from DM decay coming from these unresolved subhalos can be 
obtained by considering the flux coming from resolved subhalos, 
which is less than 1% of the flux coming from the smooth com- 
ponent. Hence, we decide to ignore the contribution of unresolved 
subhalos to the amplitude of the galactic DM-decay emission. 

Regarding the contribution to the APS from unresolved subha- 
los, we note that although subhalos just below the mass resolution 
of Aq-A-1 (~ 10 5 M Q ) might still contribute to the anisotropics, 
mainly through a Poisson-like APS, their abundance is so large 
(the subhalo mass function grows as oc M~' 9 ) that the intrinsic 
anisotropics of the gamma-ray intensity produced by them would 
be very small. Because of this, the APS at multipoles above / ~ 100 
is likely dominated by subhalos with masses above 10 5 M o (see Sec. 
I5.2l and the top panel of Fig. 8 o flAndol2009l) , allowing us to neglect 
the contribution of subhalos with lower masses. We have verified 
this is indeed the case using the analytical model of Ando ( 2009) 
(see discussion in Sec. 15. 2. 2"! and AppendixlDl. 

We do not take into account annihilation boosts due 
to fine-grained phase-space structures like streams and caus- 
tics. For a standard DM model without specific boost mech- 
anisms (e. g. Sommerfeld enhancement) these effects are sub- 
dominant ( Vogelsberger et al. 2007; White & Vogelsberger 120081 ; 
IVogelsberger et alJl2009llVogelsberger & Whitell2010h . If a mech- 
anism like the Sommerfeld enhancement is invoked, fine-grained 
streams increase significantly the main halo annihilation, but 
their contribution is typically still less than that from subhalos 
JZavalaetalfeoill) . 

Finally, in Tab. [T] we summarize the nomenclature used to 
identify the different components of DM-induced extragalactic and 
galactic emission introduced in the present section and in the pre- 
vious one. 



11 The formalism by [Kamionkowsk i et alj feplOl) overestimates the sub- 
halo abundance in the inner region of the MW-like halo, namely within 20 
kpc. To correct for this, we assume that within this radius, the spatial distri- 
bution of unresolved subhalos follows the AQ distribution, being well fitted 
by an Einasto profile with a = 0.678 and r_2 = 199 kpc. 



Name 


Description 




DJV1 halos and subhalos in Jvlo-11 catalogs 




with more than 100 particles (i.e. with 




a mass larger than M rcs = 6.89 X l0 9 M o /h) 


no-uiNivEoiviairi 


extragalactic DM (main) halos with a mass 




Knhj/nnn A/I 1 A J — fi OQ \; 1 f\8 A/f 11, 

Deiween iWnun anu ivi^q^ — o.oy x iu iviq/ii 




resolved and unresolved (sub)halos down to 




M mm . The unresolved subhalos are simulated 


EG-LOW 


following Sanchez-Conde et al. (2011) with 




k = 7x 10~ 3 




(includes EG-MSII and EG-UNRESMain) 




resolved and unresolved (sub)halos down to 




M mln . The unresolved subhalos are simulated 


EG-HIGH 


following Sanchez-Conde et al. (201 1) with 




k = 0.15 




(includes EG-MSII and EG-UNRESMain) 




smooth MW DM halo, parametrized by an 


JV1W smooth 


Einasto profile as in Navarro et al. (2008), 




and normalized to a local DM density of 




0.3 GeV/cm 3 




DM subhalos in the AQ catalogs 




with more than 100 particles (i.e. with a 




m5icc 1 Qr(Tf»r t\t an 1 "71 V 1 fl-' /W „ ^ 
llla.33 ld.lt CI Lllo.il 1 . / 1 s\ LVJ IVlQ) 




DM subhalos with a mass between M m j n 


UAL- 1 IN Kllv> (LUW ) 


and 1.71 X 10^ Mq, simulated 




following Sanchez-Conde et al. (2011) 




withyt = 7 x 10~ 3 




DM subhalos with a mass between M m j n 


GAL-UNRES (HIGH) 


and 1.71 X \0 s M o , simulated 




following Sanchez-Conde et al. (2011) 




with k = 0.2 


GAL-LOW 


MW smooth + GAL-AQ + 




GAL-UNRES (LOW) 



GAL-HIGH 



Table 1. Summary table 
the different components 



MW smooth + GAL-AQ + 
GAL-UNRES (HIGH) 

of the nomenclature used in the paper to identify 
of the DM-induced emission. 



5 ENERGY AND ANGULAR POWER SPECTRA OF THE 
DARK-MATTER-INDUCED GAMMA-RAY EMISSION 

Before showing the analysis of our simulated maps, we note that 
changing the particle physics scenario (i.e. considering a different 
value for m x and/or selecting a different annihilation/decay chan- 
nel) would require, in principle, re-running our map-making code 
for the extragalactic intensity, since the photon emission spectrum 
is redshifted along the line of sight. This is a computationally ex- 
pensive task given that one complete realization takes approxi- 
mately 50000 CPU hours. However, this is not necessary since it 
is possible, given a reference all-sky map obtained for a particular 
particle physics model, to derive the corresponding map for a dif- 
ferent model simply applying a set of re-normalization factors for 
different redshifts. Such prescription is described in detail in Ap- 
pendix [C] 
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Figure 5. Average of the extragalactic gamma-ray intensity per unit of redshift as a function of redshift at 4 GeV from DM annihilation (left panel) and DM 
decay (right panel) for \b\ > 10°. Solid black lines correspond to the contribution from resolved (sub)halos in the MS-II (EG-MSII), while the solid green lines 
include in addition the boost from unresolved main halos (EG-UNRESMain, see Section lT2l . The solid red and blue lines include all the previous components 
and the emission from unresolved subhalos down to a minimum mass M m \ n = 10~ 6 M Q according to the method described in Section 1331 for the LOW and 
HIGH case, respectively. In all cases, annihilation or decay into bottom quarks is assumed: for annihilating DM, m x = 200 GeV and (<x ;mn i>) = 3 X 10 26 cm 3 s~' , 



while for decaying DM, 



2 TeV and r = 2 X 10 s. The photon yield receives contributions from prompt emission and IC of the CMB photons. The 



dashed grey line shows the "astrophysicaf 
in Eqs.[T]and|2] 



part of the signal (with an arbitrary normalization) for the EG-MSII component, by neglecting the dN 7 /dE factor 



5.1 Analysis of the energy spectrum 

5.1.1 Extragalactic emission 

Fig.[5]shows the average DM-induced gamma-ray intensity per unit 
redshift of our simulated extragalactic maps as a function of red- 
shift (left and right panels for DM annihilation and DM decay, re- 
spectively), for an energy of 4 GeV. The average is computed over 
the whole sky except for a stri p of 10° along the g alactic plane, 
since this is the region used in lAbdo et iD {2010c) to determine 
the Fermi-LAT IGRB energy spectrum. Note that the intensity in 
each concentric shell filling up the volume of the past light cone 
is divided by the width of the particular shell in redshift space Az: 
this is roughly equivalent to computing the average of the integrand 
of Eqs. [T]and |2] over the redshift interval of each shell. The inten- 
sity from extragalactic resolved (sub)halos in the MS-II (EG-MSII) 
is shown with a solid black line. This same contribution is shown 
with a dashed grey line once the photon yield dN y /dE is removed 
from the intensity (arbitrary normalization) in Eqs.Q]and[2] leaving 
only the "astrophysical" part of the signal. In the case of annihila- 
tion, the grey line is essentially flat, with all redshifts contributing 
equally to the gamma-ray intensity (see also Fig. 1 of lAbdo et al.l 
l2010ah . Note that, in principle, the EBL attenuation should be vis- 
ible in the shape of the grey dashed line, but at 4 GeV its effect is 
negligible and the line only depends on how the DM distribution 
changes with z. In the case of decaying DM, the astrophysical part 
of the signal drops more quickly with redshift since it is propor- 
tional to the DM density (which in average grows as oc (1 + z) 3 ) 
instead of to the density squared. Once the modulation of the pho- 
ton yield dN y /dE is included, we see that the majority of the signal 
comes from low redshifts (more so for decaying DM): in order to 



contribute to the emission at 4 GeV, photons coming from higher 
redshifts need to be more energetic, and their intensity is damped 
due to a lower photon yield. For the benchmark shown in Fig. [5] 
the signal drops by a factor of ~ 3 - 5 from z = to z = 1 . 

Once the EG-MSII component is boosted up to include the 
contribution of unresolved main halos (EG-UNRESMain) with 
masses down to M m j„ = 10~ 6 M o //7, the signal increases by a factor 
of ~ 7 (~ 1.5) in the case of DM annihilation (decay). The con- 
tribution of unresolved main halos is given by integrating F mm (M) 
and FdecayCM) in Eq. [8] from M m j„ to M res . These cumulative lu- 
minosities are ultimately connected to the halo mass function and 
the single-halo luminosities L ann (M) and L dccay (M) in Eqs.|4]and|6] 
Interestingly, they combine to produce a mass-dependent contribu- 
tion that diverges towards lower masses in the case of DM annihi- 
lation (F mm (M) oc M~ 104 ), but converges in the case of DM decay 
(fctecay(A0 K AT 092 ). This is the reason why the EG-UNRESMain 
component is much larger than the resolved component in the case 
of annihilating DM, while the two remain rather similar for decay- 
ing DM. This implies that for the case of decay, the signal is es- 
sentially independent of M miu , as long as M min is low enough (see 
below). 

The total emission is obtained by summing the previous com- 
ponents and the contribution of unresolved subhalos down to M m j„. 
The LOW (red line) and HIGH (blue line) scenarios in the left panel 
bracket the uncertainty associated with the subhalo contribution, 
for a fixed value of M m j„ = 10~ 6 M o //j. We can see that unresolved 
(sub)halos boost the signal by a factor between 25 and 400 com- 
pared to the EG-MSII component. As noted before, such uncer- 
tainty is not present in the case of decaying DM, since the contri- 
bution of unresolved (sub)halos is essentially negligible. Note that 
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Figure 6. Average of the gamma-ray intensity coming from DM annihilation (left) and DM decay (right) as a function of observed energy for \b\ > 10°. Solid 
lines are for the extragalactic contribution, while dashed lines are for the galactic one. The color coding for the solid lines is the same as in Fig. \5\ while for 
the dashed lines, the green one indicates the contribution of the smooth MW halo, the black one is for resolved subhalos (GAL-AQ) and the red and blue lines 
indicate the emission from the MW smooth halo and its unresolved subhalos (GAL-UNRES) in the LOW and HIGH case respectively (only for the left panel). 
The observational data points with error bars refer to the measurement of the IGRB as given in lAbdoetaljfeoiOd) . 



the subhalo boost is smaller at high redshifts since the number of 
massive resolved main halos decreases with redshift and hence, the 
overall subhalo boost decreases as well. 

Fig. [6] shows the energy spectrum of the average amplitude 
of the extragalactic (solid lines) DM-induced gamma-ray intensity. 
We only consider an energy range between 0.5 GeV and 50 GeV, 
approxim ately the same rang e where the IGRB Fermi-LAT data are 
available JAbdo e t al. 2010c). As in Fig. [5] the average is computed 
in the region with \b\ > 10°. The left (right) panel is for annihilating 
(decaying) DM. The color-coding of the solid lines is the same as 
in Fig.|5] The full extragalactic signal, including resolved and unre- 
solved (sub)halos, is expected to lie between the solid red and blue 
lines. 

In the case of DM annihilation, the extragalactic contribution 
is dominated by unresolved (sub)halos. This prediction agrees well 
with thos e from previous wo rks. For instance, the grey band in 
Fig. 2 of lZavalaetal] d201lh can be compared with our "uncer- 
tainty" range bracketed by the red and blue line£f|. To be precise, 
the methodo l ogy implemented in the present paper and the one in 
IZavala et all < l201ll) are not identical, since the emission of unre- 
solved subhalos is accounted for in a different way. Nevertheless, 
we find that the range covered between our LOW and HIGH sub- 
halo b oosts is s imilar to those re ported in Fig. 2 of IZavala et al.l 
< l201ll) (see also lAbdo et alj2010af) . 

5.1.2 Galactic em ission 

In Fig. [6] we also show the galactic DM gamma-ray intensity, re- 
ceiving contributions from the resolved subhalos of the Aq-A-1 

12 Note, however, that although the DM particle mass and the annihilation 
chann el are the same, the annihilation cross section in Fig. 2 oflZ avala et alj 
1 201 1 ) is a factor of 5 lower than the one we use in Fig. [5] 



halo (GAL-AQ, black dashed line), the smooth MW halo (green 
dashed line), and from unresolved subhalos (down to M rain = 
10~ 6 M Q /h, red and blue dashed lines, for annihilating DM). We 
see that the emission from resolved galactic subhalos is essentially 
negligible, indeed being roughly two orders of magnitude smaller 
than the one from the smooth component (both for annihilating 
and decaying DM). The effect of unresolved subhalos is impor- 
tant only for DM annihilation and it is estimated to be between less 
than a factor of 2 (LOW, dashed red) and 10 (HIGH, dashed blue) 
times more than the smooth component. This represents an impor- 
tant difference with respect to what is found for the extragalactic 
case, where the subhalo boost can be even larger than two orders 
of magnitude. It can, however, be understood by noting that for the 
extragalactic case a given main halo and its subhalos are located 
essentially at the same distance from the observer, while for the 
galactic case, the observer is located much closer to the GC than 
to the bulk of the subhalo emission (on the outskir ts of the halo). 
This i s something that has already been noted by ISpringel et al.l 
(2008), where the subhalo boost to the smooth component of the 
Aq-A-1 halo (down to M rain = 10~ 6 M e /h) was estimated to be 
1.9, whereas for a distant observer it was 232. The value of 1.9 
is smaller than what we find for the HIGH case, even if the total 
boost of 232 for the case of a distant observer is compatible with 
our value. This is due to the slightly different radial distribution of 
the unresolved subhalos in the HIGH scenario, compared to what 
is found in lSpringel et alj {2008). 

In the case of decaying DM, the gamma-ray intensity is dom- 
inated by th e smooth compone nt (approximately compatible with 
the results of llbarra etalfeoid) . 

Comparing the total galactic and extragalactic contributions, 
we see that they are of the same order for the energy range and an- 
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nihilation/decay channel explored in Fig. This is roughly con- 
sistent with what has been reported previo usly (e.g, see Fig . 3 of 
lAbdoetalJl2010j, and also Figs. 1 and 2 of lHutsi et all201(J) . 

For the particular DM candidates explored in Fig. [6] the total 
DM-induced emission is able to account for the observed IGRB 
intensity if the HIGH subhalo boost is considered (at least in one 
energy bin). 

5.2 Analysis of the angular power spectrum of anisotropics 

We consider now the statistical properties of the anisotropies of 
our simulated maps, which is the main objective of the present pa- 
per. Two slightly different definitions of the APS will be used: () 
the so-called "intensity APS" (Q), defined from the decomposi- 
tion in spherical harmonics of the two-dimensional sky map after 
subtracting the average value of the intensity over the sky region 
considered: 



C, 



1=0 m =-l 



dE 



1 



dEj 

Z \Otm\ 



\ III ~[ 



(13) 



and ii) the so-called "fluctuation APS" (Cj? uct ), which is dimension- 
less and is obtained from the decomposition of the relative fluctua- 
tions of an all-sky map. The fluctuation APS can be obtained from 
the intensity one simply dividing by (d<t>/dE) 2 . 

The intensity APS has the advantage of being an additive 
quantity, meaning that the intensity APS of a sum of maps is the 
sum of the intensity APS of each individual component (assuming 
that the maps are uncorrected, otherwise their cross-correlations 
should also be taken into account). On the other hand, the fluctua- 
tion APS of multiple components can be summed only after mul- 
tiplying by the square of the relative emission of each component 
with respect to the total: 



i ^dE' e (a 



(d<t>'/dE) 



(d$>/dE) 



2 



(14) 



In order to compare directly the APS from our maps with the 
Fermi-LAT APS measure ment, it would be necessar y to consider 
the same target region as in lAckermann et ail d2012al) . masking out 
the known point sources and the region along the galactic plane 
(\b\ < 30°), where the contamination due to the galactic foreground 
emission is larger. In this work we only present the APS as obtained 
directly from our maps and leave the comparison to the Fermi-LAT 
APS data for future work. 

We use HEALPix to compute the APS of our template maps, 
and note that the APS is conventionally plotted once multiplied by 
t(t + l)/2n, which for larg e multipoles is pr oportional to the vari- 
ance of Aflm (see Eq. 35 of Zavala et al. 2Q10h. 



5.2.1 Extragalactic APS 

The upper panels of Fig.|7]show the fluctuation APS of our template 
maps at an observed energy of 4 GeV for the case of annihilating 
DM (left panel) and decaying DM (right panel), using the same par- 
ticle physics benchmark models used in Figs. [5] and [6] (defined in 



Sec.[2j. The color-coding is also the same as Fig. [5] solid lines indi- 
cate extragalactic components, while dashed ones stand for galactic 
ones. The minimal halo mass is assumed to be M min = 10~ 6 M o /h. 

The fluctuation APS (upper panels) illustrates clearly the dif- 
ference in the intrinsic anisotropies pattern of the different compo- 
nents which can be summarized as follow iT 4 ! 
Resolved (sub)halos in MS-II (EG-MSII): in the case of DM anni- 
hilation, the extragalactic signal from the resolved (sub)halos (solid 
black line) is less steep than a pure shot-noise power spectrum, 
characteristic of perfectly unclustered sourc es and, not surpris ingly, 
it is in agreement with the results found bv lZavala et alj d2010h (see 
the black solid line of their Fig. 12). At large multipoles, this com- 
ponent i s approximately com patible also with the top right panel of 
Fig. 2 of lAndoetalJd2007bl) . 

Unresolved subhalos of MS-II main halos: the solid yellow and pur- 
ple lines correspond to the case in which the emission of resolved 
main halos is boosted up by the contribution of unresolved subha- 
los, for the LOW and HIGH subhalo boosts, respectively. We see 
that at large angular scales, where the APS is related to the cluster- 
ing of main halos, the yellow and purple lines have a larger normal- 
ization than the black one, although their shapes are approximately 
the same. This is because subhalos give a larger boost to the most 
massive halos, which are also more clustered (biased). At interme- 
diate scales, from C = 30 to 100, the APS gets shallower reflecting 
the internal distribution of subhalos within the largest halos, which 
is considerably less peaked than their smooth density profiles. Fi- 
nally, at larger multipoles (I > 100), the emission is dominated by 
low-mass main halos and thus the yellow and violet solid lines are 
essentially on top of the solid black line. 

Unresolved main halos ( EG-UNRESMain): on the other hand, the 
solid green line indicates the case in which the contribution from 
unresolved main halos is added to the resolved component. The 
fluctuation APS of the EG-UNRESMain component alone is char- 
acterized by a lower normalization than the solid black line, since 
we assume that unresolved main halos have the same distribution of 
the least massive halos in MS-II (see Sec. l3.2t . Moreover, these are 
mainly point sources (and very numerous), thus their APS is less 
steep than the case of the EG-MSII component, being mainly sensi- 
tive to what is called the "2-ha lo term", i.e., to correlati ons between 
points in different halos (e.g. lAndo & Kom atsu 2006) . The green 
line ca n be compared with the dashed line in Fig. 12 of Zavala et all 
fcOlCh: we not e asignificant difference for £ > 40, where the APS 
in lZavala et al.1 J201fj) is closer to a pure shot-noise behaviour. This 
difference alread y appears in Fig. [2] where the APS obtained with 
the code used in lZavala et alj j2010h exhibits more power at large 
multipoles than what we find with our improved map-making code. 
We speculat e that the steep APS of the dashed line in Fig. 12 of 
IZavala et alj ( bold) is a consequence of the spurious features that 
can be seen in Fig.[2]and that we have reduced in the present work. 
Total extragalactic emission: once the unresolved subhalo boost is 
applied to halos below and above the MS-II mass resolution, we ob- 
tain the full extragalactic emission, for either the LOW (solid red 
line) or HIGH (solid blue line) subhalo cases. The contribution of 
unresolved halos (even with the subhalo boost) to the fluctuation 
APS is subdominant and the shape of the solid red and blue lines 
is exactly the same as the solid yellow and purple lines, respec- 
tively. The decrease in the normalization is due to the fact that the 



13 Notice the slightly different shapes of the energy spectra of the extra- 
galactic and galactic components due to redshifting and photon absorption 
at high energies in the case of extragalactic objects. 



14 We remind the reader that the extragalactic APS is affected by a deficit 
of power at large angular scales due to finite size of the MS-II box. 



© 0000 RAS, MNRAS 000, 000-000 



14 



Fornasa et al. 





Multipole Multipole 



Figure 7. Upper panels: Fluctuation APS of the template gamma-ray maps at an observed energy of 4 GeV for annihilating DM (left) and decaying DM 
(right). The particle physics parameters (including M m j n ) as well as the color coding are the same as those in Figs. [3] and [6] Solid (dashed) lines indicate 
the extragalactic (galactic) emission. Bottom panels: The same as the upper panels but for the intensity APS (see Eg. 113) . The upper panels give a measure 
of the relative anisotropics of the different components, whereas the bottom panels are an absolute measurement of the anisotropics and clearly show which 
components dominate the APS. The grey dashed line (with arbitrary normalization) indicates a Poissonian APS independent on multipole. 



resolved structures generate anisotropies that only contribute to a 
small fraction of the total emission (the f, factor in Eq.ll4t. 

In the lower panels of Fig.[7]we show the intensity APS, which 
allow us to estimate the absolute contribution of the different com- 
ponents. Large values of the intensity APS can be obtained from a 
particularly anisotropic component or from a very bright one. The 
angular dependence for all components is the same as in the fluctu- 
ation APS, but now, due to a very small average intensity, the EG- 
MSII component has the lowest intensity APS (black solid line), 
followed by the solid green line, corresponding to the sum of the 
EG-MSII and EG-UNRESMain components (even if the fluctua- 
tion APS is larger for the former than for the latter). Once the full 



extragalactic emission is considered (solid red and blue lines), the 
intensity APS is between a factor of 100 and 5 x 10 4 larger than the 
intensity APS of EG-MSII, depending on the subhalo boost used. 
Notice that the solid yellow and purple lines (that only include re- 
solved (sub)halos and the subhalo boost to the resolved main halos) 
have essentially the same intensity APS as the solid red and blue 
lines, which implies that the total intensity APS of the DM annihi- 
lation signal is dominated by the extragalactic unresolved subhalos 
of the massive main halos. 

In the case of DM decay (right panels), we can see that the 
fluctuation APS of the EG-MSII component (solid black line) has 
the same shape as the solid green line (which adds the contribution 
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Figure 8. Ratio of the fluctuation APS of the extragalactic maps (resolved 
(sub)halos, EG-MSII) between the case of annihilating and decaying DM 
(for the same particle physics models as in previous figures). The black line 
corresponds to the DM-induced emission up to z < 2.07 while the red line 
only accounts for the emission in the first shell (z. < 0.01). 

of EG-UNRESMain), but a higher normalization. This is because 
the signal is dominated by the massive resolved (sub)halos. We also 
see this in the case of the intensity APS (bottom right panel), where 
the contribution of low-mass halos to the intensity APS is essen- 
tially negligible (the solid green line overlaps with the solid black 
line). 

The fluctuation APS of the extragalactic maps for the case of 
DM annihilation and DM decay are very similar. This can be seen 
more clearly in Fig. [8] where we plot the ratio of the fluctuation 
APS of the EG-MSII component in the case of annihilating and de- 
caying DM. The black line corresponds to the APS of the past-light 
cone up to z < 2.07, while for the red line we only consider the 
first concentric shell (z < 0.01). The red line shows that the an- 
nihilation and decay cases are different mainly at large multipoles 
(£ > 50 - 60) where the APS is sensitive to the inner halo pro- 
file: the more extended the surface brightness profile is, the less 
steep the APS is. Thus, we expect the APS to be steeper (at large 
multipoles) for the case of annihilating DM than for decaying DM. 
However, this effect is only evident for the objects that are closer 
to us: (sub)halos that are further away appear point-like (for the 
angular resolution of the maps) and, in that case, the signal from 
annihilation and decay becomes indistinguishable, as is shown by 
the black line in Fig. [8] 

5.2.2 Galactic APS 

The dashed lines in Fig. [7J indicate our results for the APS of the 
galactic components. They can be summarized as follows: 
MW smooth halo: since the position of the observer is offset with 
respect to the GC, the DM-induced emission associated with our 
own galaxy is larger when looking towards the GC. This creates 
a large scale dipol^El that can be seen in the APS of the smooth 
component (dashed green lines in Fig. [7}, which decreases more 
rapidly for the case of decay than for annihilation since the lumi- 
nosity profile is more centrally concentrated in the latter. 
Resolved AQ subhalos (GAL-AQ): in contrast to the previous case, 
the emission from the resolved subhalos (GAL-AQ, dashed black 
lines) is much more anisotropic at larger multipoles, being rather 
similar in shape to the extragalactic one. The exact shape of the 
GAL-AQ contribution can be affected by the position of the ob- 
server relative to the local subhalo population: if a subhalo is very 
close to the observer, it would appear as a very extended source 

15 Strictly speaking the effect of having an emission peaking towards one 
particular direction does not affect only the APS at t = 1 , as a real dipole, 
but extends to much larger multipoles. 



in the sky map increasing the power at low multipoles, making the 
APS steeper. In order to quantify this effect, we constructed 100 
different sky-maps of the GAL-AQ component, randomly chang- 
ing the position of the observer in the surface of a sphere centered 
in the GC with a radius of 8.5 kpc. We find that the first and third 
quartiles of the distribution (at £ = 200) are located only a factor of 
2 below and above the median, respectively. 
Unresolved galactic subhalos (GAL-UNRES): to evaluate the con- 
tribution of unresolve d galactic su bhalos to the APS we follow the 
method presented in lAndd ( 2009), which uses analytical relations 
to calculate the APS from galactic substructures for a specified sub- 
halo distribution and luminosity function. The details of our imple- 
mentation are des cribed in App endixlDl Basing our subhalo models 
on the results of lSpringel etal] J20081) (both for a LOW and HIGH 
boost), the contribution of unresolved galactic substructures to the 
intensity APS is small: for annihilation, the contribution to the APS 
is less than ~ 10% of that from resolved subhalos, while for decay 
their contribution is at most a few percent of that from resolved 
subhalos. We therefore choose to not include this contribution to 
the APS. 

Overall, considering the galactic and extragalactic contribu- 
tions, the APS signal is clearly dominated by the smooth halo com- 
ponent in the case of DM annihilation, although the extragalactic 
emission could be important at very large multipoles (£ > 300) if 
subhalos give a large boost. On the contrary, for DM decay, the ex- 
tragalactic emission dominates already from £ > 20 and it is only 
at the very large scales that the anisotropy of the smooth halo dom- 
inates the sign al. However, if a mask is introduced along the galac- 
tic plane (as in lAckermann et alj2012ah . we expect that the balance 
between galactic and extragalactic components will change, reduc- 
ing significantly the impact of all the components characterized by 
a large emission around the GC (see Sec.[6]l. 



6 DISCUSSION 

6.1 Redshift dependence of the extragalactic APS 

In Fig. [9] we divide the extragalactic gamma-ray emission in red- 
shift bins (each bin including 4 MS-II snapshots) and compute the 
fluctuation APS for the EG-LOW component (i.e. the total emis- 
sion, including all (sub)halos, down to M min = 10~ 6 M G //i and with 
a LOW subhalo boost) in each bin. The APS is computed at an 
energy of 4 GeV. We can see that for both, DM annihilation (left 
panel) and DM decay (right panel), the lower redshifts are charac- 
terized by a larger anisotropy. This is due to the fact that the volume 
of the past light cone grows with redshift, as well as the number of 
gamma-ray emitting (sub)halos. Thus, the first snapshots are those 
characterized by the lowest number of (sub)halos and are more 
affected by their discrete distribution. Moreover, the clustering of 
DM (sub)halos is larger at lower redshifts. The peaks that move to- 
wards higher multipoles with increasing redshift are a remnant of 
the spurious effect related to the periodicity of the MS-II box dis- 
cussed in Sec. 13.11 For a particular redshift, the peaks indicate the 
angular size of the MS-II box at that redshift (what we called I* in 
Sec. l3.lt : multipoles smaller than I* are affected by a loss of power 
due to the missing modes at wavelenghts larger than the MS-II box, 
and therefore we cannot trust our predictions below £* . This fact is, 
however, not relevant for a comparison with the Fermi-LAT APS 
data, since we are mainly interested in the multipole range between 
£ = 155 and 500, where the extragalactic APS is dominated by the 
first redshifts, for which £* < 20 - 30. 
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Figure 9. Fluctuation APS for the extragalactic gamma-ray intensity coming from DM annihilation (left) and DM decay (right) for different redshift bins. The 
APS is computed at an energy of 4 GeV and for the LOW subhalo boost (in the case of annihilation) with M m j n = 10~ 6 M e /h. The fo-model is assumed. 



6.2 Energy dependence of the APS 

The extragalactic fluctuation APS incre ases with increasing en- 
ergy, a fact already pointed out in the past dAndo & Komatsul2006l : 
IZavala etal . 2010; Iba rra et alj2010h and related to the redshift de- 
pendence discussed in the previous section: following Eq. [14] the 
total fluctuation APS at a particular energy E y can be written as 
the sum 2, f^(E y )Cf" a over the fluctuation APS of each concentric 
shell Cf uct normalized by the square of the relative emission in the 
;-th shell with respect to the total. Since individual shells are thin 
in redshift space, each single Cf uct does not depend on energy, and 
thus changing the energy only has the effect of modifying the /; 
factors that determine the balance among the APS of the different 
shells. These f, factors depend on the annihilation/decay channel 
selected for the particular DM candidate, as well as on how much 
the DM density changes with z within a particular shell (see Fig.[5]l. 
For high energies, the shells that contribute the most to the signal, 
i.e. those with the largest fi factors, are the first shells, which are 
characterized by the largest APS. Thus, the total fluctuation APS 
increases as energy increases. 

Note, however, that in the cases where the fluctuation APS is 
dominated by the galactic emission, the fluctuation anisotropy will 
not change with energy, neither in normalization nor in shape. 



6.3 Inner density profile of DM (sub)halos 

When dealing with the extragalactic emission, we have assumed 
that (sub)halos have a smooth NFW density profile, which is char- 
acterized by a slope that tends asymptoticall y to -1 for small radii, 
and lies between the st eeper Moore (Moore et al. 1999) and cored 
Burkert (Burkert 1996) profiles. Current high resolution W-body 
simulations have demonstrated that the Einasto profile (Eq.ll2t pro- 
duces an even better fit than NFW. The slope of the Einasto profile 
decreases as a power-law as the distance from the center decreases, 
and it is shallower than NFW at small radii. 

It is also important to note that the process of galaxy forma- 
tion within DM halos has an impact on the DM distribution in the 
central regions where it is believed that the halo is adiabatically 
contracted resulting in a more concentrated DM distribution (e.g. 
iMo et alJ[l998l ; lGnedin etal]|2004lAnn et al.ll2007h . However, re- 



cent hydrodynamical simulations with strong supernovae feedback 
claim that including the effect of baryons can actually result in 
the development of a centra l DM core in interme diate mass halos 
jPontzen & Governatdl201ll ; lMaccio' et alfeoill) . 

For the extragalactic emission, the uncertainty on the inner 
DM density profile has very limited effect on the APS since only a 
small fraction of (sub)halos cover more than one pixel in our maps, 
and also because, even if the object is characterized by extended 
emission, the difference between a cuspy or cored profile is only 
noticeable at very small projected radii. On the other hand, we do 
expect a change in the total intensity of the DM-induced emission: 
for instance, if an Einasto profile is used ins tead of a NFW, th e 
annihilation rate per halo will increase by 50% dZavala et alj20ld) . 
Considering the extreme cases of a Burkert and a Moore profile, the 
differ ence is roughly an order of magnitude ( tProfumo & Jeltem3 
2009). For decaying DM the total luminosity of a halo is directly 
proportional to its mass, independently of the DM profile assumed. 

For the galactic emission, the reasoning above applies to the 
resolved subhalos. On the other hand, assuming a different profile 
for the smooth halo may have a stronger impact on the APS, since 
this represents the largest contribution (at least at low multipoles, 
and particularly for the case of annihilating DM). The effect of as- 
suming a NFW3 rather than an Einasto profile is evident at low 
multipoles (with the APS of the former being smaller than the APS 
of the latter), but the difference becomes smaller at larger multi- 
poles. This can be explained by noting that the emission towards 
the GC is lar ger with respect to the anti-center in the case of an 
Einasto profilq^J, resulting in a more anisotropic APS. The differ- 
ences are less evident for the case of decaying DM. 

Finally, it is important to remember that any uncertainty in the 
inner MW density profile will be reduced if the region around the 
GC is masked. For instance, for \b\ > 30°, th e different reasonabl e 
DM profiles are practically indistinguishable teertone et al.11200^) . 

16 Taken from fPrada et alj|2004l and normalized to the same local density 
than the Einasto profile introduced in Sec. 14. II 

17 If the two profiles are normalized to the same density at 8.5 kpc, the 
NFW will have a larger intensity within ~ 10 pc, but in the region between 
~ 1 kpc and 10 kpc from the GC (where the majority of the emission comes 
from) an Einasto profile is characterized by a larger density. 



© 0000 RAS, MNRAS 000, 000-000 



Characterization of Dark-Matter-induced anisotropics in the diffuse gamma-ray background 17 










i i _ 




M m 


n =10" 12 luL/h 




: 


M m 


„=1<y 6 NUh 






M m 








M m 


n =10 3 W^/h 






M m 


„=10 6 wjh 






M m 


n =6.89 10° MJh ^ 




M m 


n =10"W^/h 





































10 



10 



Multipole 



10 d 



W 
Q_ 
< 

o 
3 



o 1 



o 
g 



10" 



^~h^V»^A^^ 


M m 

- M m 

M m 

M m 

M m 

M m 

M m 

I I 


n =10" 12 WL/h 
n =10' 6 WL/h I 
„=1 Mj/h 
n =10 3 Wlj/h 
n =10 e WL/h 
„=6.89 10 s tJ^h - 
n =10 12 ^h 



10 



Multipole 



10" 



Figure 11. Ratio of the total fluctuation APS (galactic and extragalactic) for different values of M mm with respect to the reference case M m 



10- 6 M o //i. 



The APS is computed at 4 GeV. The left panel refers to the case of annihilating DM with the HIGH subhalo boost, while the right panel is for decaying DM. 
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Figure 10. Total gamma-ray intensity from DM annihilation (solid lines) 
and DM decay (dashed lines) at 4 GeV as a function of the minimal halo 
mass M m j„. The LOW and HIGH subhalo boosts (see Secs. l3.3 land l4.21 are 
shown with black (red) lines, respectively. The emission has been computed 
only for the values of M ra j n indicated by the full dots, while the lines are 
obtained by interpolation. 



6.4 The minimum self -bound halo mass M mill 

The particle nature of DM determines the small-scale cutoff in 
the matter power spectrum of density fluctuations, and hence, the 
value of M min . For neutralinos, the most common WIMP DM can- 
dida tes, typical values for M min go from l0~ l[ M e /h to 10~ 3 M o /7i 
(e.g. lBringmann| r2009). Although this range can be considered as 
a reference for all WIMP candidates, a particular scenario might 
lie outside this range. In order to investigate the impact of differ- 
ent values of M mm in our predictions, we generate template maps 
for M m i n equal to 10~ 12 M o //j and 1 M Q /h. We also consider a few 
larger values (namely M min = 10 3 , 10 6 , 6.89 x 10 8 and 10 I2 M o //z) 
that, although clearly far above the expected mass range for WIMP 
models, are discussed in order to understand how halos of different 
masses contribute to the gamma-ray intensity and APS. 

In terms of the mean gamma-ray intensity, we can see the 
impact of changing M m ;„ in Fig. 1101 solid (dashed) lines for the 



case of DM annihilation (decay), while black and red lines refer 
to the LOW and HIGH subhalo boosts, respectively. For annihi- 
lation, the mean flux decreases only by a factor of ~ 5 between 
M mm = lQ- ,2 M Q /h and 1 MJh, for the LOW case, while the dif- 
ference is one order of magnitude for the HIGH case. For even 
higher values of M mm , the intensity stays essentially constant for 
the LOW case, while it decreases further for the HIGH case un- 
til reaching a plateau. In both cases, the point where the intensity 
reaches an approximately constant value marks the region where 
the total intensity passes from being dominated by the extragalac- 
tic component (for low M ra j n ) to being dominated by the emission 
in the MW (for larger M rain ). The transition happens at the larger 



values of M mm for the HIGH subhalo boost, because (see Sec. 15. lb 
increasing the subhalo abundance produces more significant effects 
for the extragalactic emission than for the galactic one. 

In the case of DM decay (as discussed before), the bulk of 
the emission is dominated by large mass halos and, in particular, is 
already accounted for in (sub)halos with masses larger than 10 s - 
l0 9 M o /h, with smaller objects contributing only marginally. 

The effect of M mm on the total APS is shown in Fig. [TT] the 
two panels indicate the ratio of the fluctuation APS at 4 GeV for 
7 values of M min with respect to the reference case of = 
l0~ 6 M o /h. The panel on the left shows the case of an annihilat- 
ing DM candidate with a HIGH subhalo boost: looking at Fig. 
[7J for M mm = 10~ 6 M o //i, the total intensity APS is dominated 
by the MW smooth halo, while the contribution of extragalactic 
(sub)halos plays a role only at large multipoles. Now, going from 
M m i n = 10~ 6 to 10~ l2 M o //i does not have a strong impact on the 
galactic component but it makes the total extragalactic emission 
increase by a factor of a few (see Fig. 1 10t. The net effect, follow- 
ing Eq.[l4] is that the total fluctuation APS decreases because less 
intensity is associated with the component that dominates the in- 
tensity APS (i.e. the galactic one). This is also the reason why the 
total fluctuation APS increases from M mi „ = lO" 6 to 1, 10 3 and 
10 6 M o //z. When the total emission starts to be dominated by the 
MW smooth halo (i.e. above approximately ~ 10 6 M o //j), there is 
essentially no change to the APS due to variations in Mmm. 

The same features appears in the case of a LOW subhalo boost 
(the figure is not present), even if this case is characterized by a 
smaller relative difference (all the lines are within one order of 
magnitude), and, since the emission of the DM smooth halo starts 
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Figure 12. Energy spectrum of the average gamma-ray intensity from DM 
annihilation (color lines) or decay (black line) from extragalactic and galac- 
tic (sub)halos. The blue and red lines correspond to the LOW and HIGH 
subhalo boosts, respectively, so that the filled grey area between them cor- 
responds to the uncertainty due to the subhalo boost, for a fixed value of 
Mmin, The red (blue) shaded area around the red (blue) solid line indicates 
the uncertainty in changing the value of M m j„ from 1CT 12 to 1 Mg/h, for the 
LOW (HIGH) scenario boost. The solid black line shows the prediction for 
a decaying DM candidate and the black shaded area (appearing as a thick- 
ening of the solid black like) indicates the uncertainty in changing M mm 
from 1CT 12 to 1 M G /h. The observational d ata points with error bars refer 
to the measurement of the IGRB as given in lAbdo et alj feOlOch . Only the 
emission with \b\ > 10° is considered. The DM candidates are described in 
Sec.E] 

to dominate already at M mm = 1 M Q /h, the APS does not change 
for M mm larger than that value. 

The right panel in Fig.QT|is f° r decaying DM: the different 
lines follow the same behaviour as for annihilating DM but the ef- 
fect of changing M rain is highly reduced. The only important devi- 
ation is for the largest value of M mm : at low multipoles the APS is 
still dominated by the smooth MW halo and, thus, we expect only 
a different normalization. However, when the extragalactic compo- 
nent becomes relevant, the dashed red line decreases because the 
extragalactic fluctuation APS for M m j„ = l0 [2 M o /h is smaller than 
the case at M min = 10~ 6 M G /7i since it is sensitive, at these multi- 
poles, to the inner DM profile of the largest objects. 



6.5 Theoretical uncertainty bands 

In the current section we summarize our predictions for the energy 
and angular power spectra of the DM contribution to the IGRB 
emission. We also present "theoretical error bands" that bracket the 
uncertainties discussed in the previous sections. These predictions 
are given only for a fixed particle physics scenario (the fo-model, 
see Sec. |2j, while the analysis of different DM candidates, (i.e., 
changing m x , the annihilation cross section, decay life time and an- 
nihilation/decay channels), will be discussed in a follow-up paper. 

The energy spectrum of the DM-induced signal (averaged over 
the region with \b\ > 10°) is shown in Fig. [12] The grey area be- 



tween the red (LOW subhalo boost) and blue line (HIGH subhalo 
boost) spans approximately a factor of 50 and quantifies the uncer- 
tainty associated with the unknown subhalo boost, for a fixed value 
of Mmm = 10~ 6 M o //j. The additional red and blue shaded areas in- 
dicate the uncertainties introduced by changing the value of Mmm 
between 10~ l2 M o /h and l M Q /h. For the case of decaying DM, 
our predictions are completely determined by massive (sub)halos 
so the theoretical uncertainties are much smal ler than for t he case 
of DM annihilation. The Fermi-LAT data from lAbdo e7aT] ( l2010ch 
are also plotted with error bars. 

Fig. 1 1 31 summarizes our predictions for the DM-induced APS 
(intensity APS in the left panel and fluctuation APS in the right 
panel). Contrary to the plots presented in the previous sections, the 
APS is now computed after having integrated the gamma-ray emis- 
sion between 2 and 5 GeV. Moreover, the APS has been averaged in 
bins of At = 50 starting from t = 5, and we introduce a mask cov- 
ering the region with \b\ < 30°. We approximately correct for the 
effect of the mask by dividing the raw APS by the fraction of th e 
sky / sky left unmasked, as it was done in lAckermann et al] J2012al) . 
All of this is for comparison purposes with the Fermi-LAT APS 
data in the same energy bin, taken from lAckermann et al. 2012^^1 
The inclusion of the mask has strong effects both on the average 
emission of the smooth MW halo and on its APS since we are 
masking the region where the signal peaks. On the other hand, it 
has a limited effect on the extragalactic emission. After masking, 
the total intensity APS for annihilating DM is dominated by the re- 
solved galactic subhalos in the case of the LOW subhalo boost and 
by the extragalactic unresolved (sub)halos for the HIGH subhalo 
boost, i.e., contrary to what is shown in Fig. [7] the smooth MW 
halo only represent a subdominant contribution. For decaying DM, 
all these three components (extragalactic emission, resolved galac- 
tic subhalos and the smooth MW halo) have a comparable intensity 
APS. 

In Fig. [T3] the red and blue lines indicate our predictions for 
an annihilating DM candidate in the LOW and HIGH scenario, re- 
spectively. Thus, the grey area indicates the uncertainty associated 
with the unknown subhalo boost. If we had plotted only the extra- 
galactic intensity APS in the left panel, the LOW case would have 
been a factor 500 below the line for the HIGH case (as in Fig. [7j. 
However, the resolved galactic subhalos increase the intensity APS 
for the LOW case, while having a less important role for the HIGH 
case. Thus, the red and blue lines are only one order of magni- 
tude away from each other. Moreover, the uncertainty due to M min 
is completely negligible in the LOW case since the APS is deter- 
mined by the galactic resolved subhalos, and thus is not sensitive 
to changes in M mm . The same is true for the case of decaying DM 
(black line), whose APS is determined by massive (sub)halos. 

The right panel of Fig. Q~3] shows the fluctuation APS: the red 
line, corresponding to the LOW subhalo boost is now above the 
blue line, relative to the HIGH subhalo boost. This is because the 
galactic subhalos (the component that dominates the total APS in 
the former case) are associated with a larger intrinsic anisotropy 
than the extragalactic (sub)halos, which dominate the APS in the 
latter case. 

We conclude this section with a comment on the comparison 
between our predictions for the DM-induced APS with the Fermi- 
LAT data shown in Fig.[T3] Although a rigorous comparison is left 
for future work we can already see that the fluctuation APS from 



We do not mask th e point sources in the 1-year catalog, as in 
lAckermann et alj J2012at) , so that our f^ y is 0.5. 
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DM annihilation is of the same order, and has a similar shape, as 
the data (at least in the case of the LOW subhalo boost). On the 
contrary, for a decaying DM candidate, the predictions are not com- 
patible with a flat APS and they are also characterized by a normal- 
ization which is too low. 

Nevertheless, even if the annihilating DM candidate we used 
here is able to reproduce the same level and shape of the fluctu- 
ation APS inferred from the data, it still does not represent a vi- 
able interpretation, since such a candidate is characterized by a very 
low intensity APS (left panel). Improvements in the analysis of the 
gamma-ray data are still possible both from the experimental side 
(e.g. increasing the statistics, especially at high energies), and from 
the theoretical side (e.g. one can think of selecting, for each DM 
candidate, the energy bin that maximizes the DM-induced intensity 
APS). 



7 SUMMARY AND CONCLUSIONS 

In the present paper we generated all-sky gamma-ray maps from 
the annihilation/decay of DM in extragalactic (sub)halos and in the 
halo and subhalos of the MW. Apart from the prompt gamma-ray 
emission, we also considered emission due to the IC scattering 
of e + /e~ produced in the annihilation or decay with CMB pho- 
tons and, for the smooth MW halo, additional contributions from 
starlight (either directly or re-scattered by dust) and the so-called 
"hadronic emission" (see AppendiceslAlandlBl are also considered. 

The DM distribution was modeled using state-of-the-art N- 
body simulations: Millennium-II for extragalactic (sub)halos and 
Aquarius (Aq-A-1) for the galactic halo and its subhalos. To com- 
pute the e xtragalactic emissio n, we improved the algorithm de- 
scribed in Zaval a et alj d201Ct) and simulated the past light cone 
up to z = 2. The MS-II allows us to account for the emission 
of structures with a mass larger than M tes ~ l0 9 M o /h. We then 
considered the intensity from unresolved (sub)halos down to a 
minimum self-bound mass M min , by a hybrid method that com- 
bines an extrapolation of the behaviour of the least massive re- 
so lved halos in MS-II with the subhalo bo ost model introduced 
in Kamionkowski & Koushiappas (2008) and Kamionkowski et al. 
d2010h and refined in ISanchez-Conde et alj j201lh . On the other 
hand, the galactic emission was modeled assuming that the smooth 
halo of the MW is given by an Einasto profile, renormalized to a 
value of 0.3 GeV/cm 3 for the local DM density. Resolved galactic 
subhalos are taken directly from the Aquarius simulation (down to 
a mass of w 10 5 M o ), while the contribution of unresolved galactic 
subhalos is estimated by means of the same procedure used for the 
extragalactic emission. 

The template maps of the DM-induced emission were then 
used to derive the energy spectrum of the different components 
from 0.5 GeV and 50 GeV (see Fig.(6](. The main goal of the paper 
is the characterization of the anisotropies of the DM-induced emis- 
sion, which was done in Sec. 15. 21 where we computed the APS of 
the different components up to C = 500 (see Fig. [7}, which is the 
range covered by the recent Fermi-LAT analysis of th e APS of the 
diffuse gamma-ray emission dAckermann et al.ll2012al) . 

We also discuss the possible effects of modifying some of the 
assumptions in our modeling of the DM distribution. Most notably, 
we consider two different scenarios with a small and a large subhalo 
contribution (referred to as LOW and HIGH throughout the text). 
Additionally, we study how the energy spectrum and APS depend 
on the value of the minimal self-bound halo mass M m j„. A discus- 
sion on the effects of using different DM halo profiles is also given. 



Quantifying the impact of these uncertainties helps us understand- 
ing which are the ones that primarily affect the APS, as well as to 
associate a "theoretical uncertainty band" to our predictions. 
The main results of our study are: 

• An improvement of the procedure used in lZavala et alj J2010t) 
to compute the extragalactic DM-induced intensity introducing in- 
dependent rotations for each of the replicas of the simulation box. 
This notably reduces spurious features in the APS of the simulated 
maps due to residual correlations introduced by the periodicity of 
the MS-II box. 

• For annihilating DM, the total extragalactic emission (once all 
(sub)halos down to M mm = 10~ 6 M Q /h are considered) is a factor 
of 20 (500) larger than the emission produced in the (sub)halos re- 
solved by the MS-II simulation. On the other hand, the extragalactic 
emission for decaying DM is dominated by the structures resolved 
in the simulation, with a total intensity that only increases by a fac- 
tor of 2 once unresolved objects are taken into account. 

• The effect of including unresolved subhalos is less important 
for the galactic component, since these are mainly located in the 
outskirts of the MW halo, far from the observer, contrary to the 
nearby GC that produces a significant contribution to the signal. 
Our prediction for the total galactic intensity (down to M m j„ = 
1O~ 6 M //?) is between a factor of 2 and 10 times larger than the 
emission of the smooth MW (for annihilating DM). The contribu- 
tion of unresolved subhalos is negligible in the case of DM decay. 

• The extragalactic intensity APS in the case of annihilating DM 
is dominated by unresolved (sub)halos. The intensity APS of the to- 
tal emission is between 100 and 5 x 10 4 times larger than if only 
the resolved MS-II (sub)halos are considered, even though its fluc- 
tuation APS is lower than the fluctuation APS of the resolved com- 
ponent. In the case of the galactic substructures, the intensity APS 
is dominated by the resolved subhalos (which have the largest in- 
trinsic anisotropies of all components) in the Aquarius halo (down 
to ~ 10 5 M o ), while unresolved subhalos are not expected to con- 
tribute. The total intensity APS is dominated by the smooth DM 
halo of the MW, at least for low multipoles, while above I = 300, 
the extragalactic contribution becomes important (if the HIGH sub- 
halo boost is assumed). 

• The case of decaying DM is quite different: the APS of the 
smooth MW halo decreases more rapidly, so that the total inten- 
sity APS is dominated by extragalactic halos around I = 20 - 30. 
Galactic subhalos, on the other hand, are characterized by large 
anisotropies but their low intensity forces them to play only a mi- 
nor role in the total intensity APS. 

• Both for annihilating and decaying DM, the total intensity 
APS depends mainly on structures in the local Universe, with ob- 
jects located at z > 0.26 contributing to less than 10% of the total 
signal. 

• Changing the value of M mj „ from 1 to 10~ 12 M o //j has a very 
small effect for decaying DM, while our predictions can change 
dramatically for annihilating DM, especially for a HIGH subhalo 
boost: the left panel of Fig.|13|shows that an uncertainty of almost 
two orders of magnitude is associated with the total intensity APS 
in this case. 

In a future work the DM template maps produced here will 
be used to derive constraints on the particle physics nature of DM 
from a comparison with the Fermi-LAT data. In Fig. [13] we made 
a first comparison for a particular DM candidate used in this work 
and find that even if the DM-induced fluctuation APS is of the same 
order of the Fermi-LAT data (for DM annihilation), this particular 
DM candidate is not able to account for the bulk of the signal de- 
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Figure 13. Total intensity APS of the gamma-ray emission from DM annihilation (color lines) or decay (black line) in extragalactic and galactic (sub)halos. 
The blue and red lines correspond to the LOW and HIGH subhalo boosts, respectively, so that the filled grey area between them corresponds to the uncertainty 
due to the subhalo boost, for a fixed value of M mm . The red (blue) shaded area around the red (blue) solid line indicates the uncertainty in changing the value 
of M min from lfT 12 to 1 M B /h, for the LOW (HIGH) case. The solid black line shows the prediction for a decaying DM candidate and the small black shaded 
area, appearing as a thickening of the solid black line, indicating the uncertainty in changing M m j n from 10 -12 to 1 M e /h. The APS is measured in the energy 
bin between 2 to 5 GeV. The observational data points with error bars refer to the measurement of the APS as given in Ackermann et al. 1 2012a). A region of 
30° around the galactic plane has been masked and the APS has been binned with a binsize of At = 50. The DM candidates are described in Sec. [2] 



tected by Fermi-LAT since its intensity APS is too low. A more rig- 
orous comparison (coupled with a scan over a reasonable set of DM 
models and using a broader energy range) is still required in order 
to derive more conclusive statements. Based on the energy spec- 
tra of the DM candidates considered here relative to the measured 
IGRB (see Fig. [6}, the APS of the 2-5 GeV energy band shown in 
Fig. [T3]is likely not the optimal choice for setting constraints, but 
it should be considered as an example for the comparison between 
the Fermi-LAT data and our predictions. 

It is also important to note that the majority of the IGRB emis- 
sion is expected to be produced by standard astrophysical unre- 
solved sources, such as blazars, star-forming galaxies and pulsars. 
Thus, a complete study of the IGRB emission can only be per- 
formed with a model that also includes these contributions. In this 
case, the so-called "energy anisotropy spectrum", i.e. the fluctua- 
tion APS at a fixed multipole but as a function of the energy, is a 
particularly useful observable since it has been shown that mod- 
ulations in the energy anisotropy spectrum may mark transitions 
between regimes where different class es of sources are responsi- 
ble fo r the bulk of the IGRB intensity dSiegal-Gaskins & Pavlidovj 
l2009h . 

In any case, the study of the IGRB energy spectrum and of its 
anisotropics are not the only tools on e can resort to f or the study 
of the IGRB nature. For instance, in IXiaet alj JioTlh the authors 
compute the cross-correlation of the Fermi-LAT data with the an- 
gular distribution of objects detected in different galaxy surveys. 
Assuming that these objects represent the detected counterparts of 
unresolved astrophysical sources contributing to the IGRB, they 
used the cross-correlation me asurement to put cons t raints on the 
IGRB com position. Moreover. iDodelson etall d2009h . [Baxter et all 
J2010h and lMalvshev & Hoed d20 1 lh showed that the analysis of 



the probability distribution of the photon counts can be used ef- 
ficiently to distinguish a DM signal from a cumulative emission 
of astrophysical sources in the IGRB data. In principle, the maps 
produced in the presen t paper represen t uni que tools to extend the 
techniques exploited in IXia et al.1 f20Tlh and lDodelson et alj (2009) 
by including a possible DM contribution. 
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APPENDIX A: INVERSE COMPTON EMISSION 

The secondary IC e missio n has been described in detail in 
iBlumenthal & Gould ( fl97(jh . where the authors also provide use- 
ful formula; to reproduce their calculation. This process consists 
in a transfer of momentum from a high energy cosmic ray (CR) 
electron or positron to a low energy ph oton of the ISRF. 

A model for the ISRF provided by Moskale nko et alj {2006) is 
publicly available on the GALPROP webpagj*^. In order to com- 
pute the IC emission using semi-analytical methods, it is conve- 
nient to fit the GA LPROP model of the I SRF as a sum of five black- 
body spectra (e.g. lDelahave et al]|2010h . One of these is the CMB, 
while the others come from a fit to the model and have less physi- 
cal meaning, although they derive from dust and stellar emissions. 
In this procedure, it is necessary to assume a homogeneous ISRF 
which might impact on the morphology of the resulting gamma- 
ray emission, although it should be quite moderate since variations 
of the ISRF affect both the e + /e~ spatial density and gamma-ray 
emissivity in opposite directions. 

Apart from the ISRF, one also needs to know the e + je~ dis- 
tribution and propagation in the galaxy in order to compute the IC 
emission. These processes are governed by the following diffusion- 
loss equation (neglecting convection and re-acceleration effects): 

- V ■ (D(E, x)V/) - —(b(E)f) = Q(E, x), (Al) 
oE 

where f(E, x) is the e + /e~ number density per unit of energy at 
the point x, D(E, x) is the diffusion coefficient while b(E) describes 
the energy losses (due to synchrotron and IC emissions). Finally 
Q(E, x) indicates the source term which in our case is DM annihi- 
lation/decay. 

Eq. I Al I governs diffusion inside a so-called diffusion zone, 
outside of which electrons and positrons are not confined by mag- 
netic fields and escape from the galaxy. The coefficients defining 
the different terms in Eo. lAll are constrained by the available ob- 
servational data (mainly the boron-to-carbo n CR ratio), but im por- 
tant uncertainties are still present (see, e.g. iDonato et al.ll2004l and 
their definition of the MIN/MED/MA X scenarios). We use h ere the 
semi-analytical methods described in Dela have et aT] d2008l) which 
take into account the full expression of the energy losses in the 
Klein-Nishina regime; 

As explained in Boehr rTet al.l ( feOlOh . the morphology of the 
galactic IC emission created by the e + je~ produced by DM anni- 
hilation/decay is very sensitive to the choice of the CR propaga- 
tion parameters and hence, the results should be taken with cau- 
tion. Here we use the same propagation model parameters as for 
the protons and anti-protons related to the hadronic emission (see 
Appendix [Bj and we assume the MED scenario mentioned above. 
The uncertainty in the resulting gamma-ray intensity can be quite 
large, depending on the arrival direction, and the results can also 
change with different e + je~ propagation models. Nevertheless, we 
neglect this source of uncertainty noting that IC emission is rele- 
vant only for a fraction of the energy range considered here and 

19 http://galprop.stanford.edu/ 



only for massive DM candidates (see Fig. QJ: in the case of the b- 
model, with a mass of 200 GeV, the IC emission is located almost 
2 orders of magnitude below the prompt emission and is dominated 
by interactions with the ultraviolet component of the ISRF. In the 
case of a decaying DM particle, though the mass is higher, the sig- 
nal gets stronger because it is not concentrated around the GC and 
the average over the whole sky is larger. Moreover, for the case of 
decaying DM, the signal is proportional to the inverse of the DM 
mass, whereas in the annihilating case it is inversely proportional 
to its square. For the T-model the same difference appears between 
annihilation and decay and it is even stronger since the masses are 
the same for both cases. Moreover, since in this case the prompt 
emission is much lower than for the i-model, IC and prompt emis- 
sion become of comparable importance, especially below 10 GeV. 

The DM-induced IC emission is implemented in a differ- 
ent way for the different components that constitute the emis- 
sio n: for the smooth halo o f the MW, the complete ISRF given 
by iMoskalenko e t al. (2006) is used, solving Eq. I Al I and consid- 
ering the propagation of e + je" produced by DM annihilation/decay 
before they interact with the ISRF. On the other hand, for the extra- 
galactic (sub)halos and for the galactic subhalos (both resolved and 
unresolved), we only consider IC scattering with the CMB photons, 
and ignore the propagation of e + je~ (spatial diffusion is ignored 
in Eq. IA1I and only IC energy losses are considered). In princi- 
ple, the secondary IC emission from massive halos (and some of 
the most massive subhalos) may be more realistically described if 
a full propagation model that includes the effect of baryons and 
secondary emission contri buted by starlight and infrared light is 
applied instead (e.g., see IColafrancesco et al. I l2006l . l2007l for the 
case of the Coma galaxy cluster and the Draco dwarf spheroidal). 
However, the contribution of extragalactic structures and galactic 
subhalos is dominated by low-mass objects where star formation 
is highly suppressed and thus are expected to have a rather small 
stellar component or be devoid of stars. 

Fin ally, we note that when w e use the complete ISRF pro- 
vided in IMoskalenko eTa l. (2006), the template maps for the IC 
emission (Cartesian maps with 90x180 pixels) have a poorer reso- 
lution than the maps obtained from the prompt emission (HEALPix 
Maps with N_side=512) due to the substantial numerical effort re- 
quired in solving Eq. I A 1 1 For our purposes it is enough to re-bin 
the Cartesian maps into a HEALPix pixelization. 



APPENDIX B: HADRONIC EMISSION 

The hadronic emission is the mechanism that contributes the least 
to our signal but it has a different spatial morphology with respect 
to the others considered in this work. We only account for it in the 
case of the smooth MW halo. It comes from the interaction of CR 
protons and anti-protons with interstellar gas. To compute such a 
component one needs to derive the pip intensity everywhere in the 
diffusion halo. To d o so, we follow the semi-analytical method of 
iBarrau et al.l d2002l) using the p ropagation parameters of the MED 
scenario in lDonato et "aT] fe004l) which gives a good fit to the boron- 
to-carbon observational data. Once the pip distribution has been 
obtained, it is con volved with the gam ma-ray production cross sec- 
tiorF°l (taken from Hua ng et"al] 2007) an d the interstellar gas^ distri- 
bution (taken from lPohl et al.l2008h . See lDelahave et alj I2OI1I) for 
more details on the computation. 

20 Note that, following iF.Wj 19671 we consider here that protons and anti- 
protons have the same cross-section. 
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Contrary to the IC emission, this component is less dependent 
on the choice of propagation parameters due to the fact that protons 
propagate much further than electrons and tend to smooth out all 
small scale effects. Moreover, the hadronic emission naturally fol- 
lows the interstellar gas distribution. A source of uncertainty, which 
we neglect, may come from the presence of DM substructures near 
the galactic disk, which may alter locally the gamma-ray intensity. 

We note that the angular resolution of the hadronic com- 
ponent is mainly limited by the resolution of the gas maps: 
0.5 deg x0.5 deg Cartesian maps. As for the case of the IC emis- 
sion, these maps are transformed into HEALPix maps. 



APPENDIX C: MAP MAKING FOR A GENERIC 
PARTICLE PHYSICS MODEL 

In this section we describe the implementation of an approximate 
method that can be used to obtain a full-sky map of the DM-induced 
extragalactic emission for any particle physics model, given a refer- 
ence map obtained for a specific model. Thanks to this method we 
only need to run once our map-making code, saving computation 
time. 

For the purposes of computing the DM-induced emission, 
each particle physics model is defined by the mass of the DM candi- 
date, its annihilation cross section (<x ann v) (or decay lifetime t), and 
the branching fractions for different annihilation (or decay) chan- 
nels with the corresponding photon yields, 6, and dNj/dE. Unless 
the model has a velocity-dependent cross section (a case we do not 
explore here), (cr ann v) is constant, as it is r in the case of decay, and, 
therefore, it is just a multiplicative factor in Eqs.Q]and[2] The pho- 
ton yield, however, depends on redshift. In the case of the galactic 
halo emission, the redshift variation across the DM sources is neg- 
ligible, but for the extragalactic component, the additional integra- 
tion over redshift links together the astrophysical and the particle 
physics factors in Eqs.Q]and[2] 

We benefit from the fact that the past-light cone in our sim- 
ulated maps is divided in concentric shells with a small redshift 
width Az (see Sec. 13. 11 1. It is then always possible to find a partic- 
ular redshift value contained within Az (called z re f) so that we can 
take the factor dN y /dE outside the integral in Eq. Q]and write the 
intensity coming from that shell aj 2 '! 

di {E ^ M) = -2> m >< 

X i 

J dA(z)p 2 (Mz), y f')e- T ^ (z ^ ) . (CI) 

In principle, each line of sight (pixel) in the sky map will have a 
different value of z rc f (for the same shell) since the integrand in the 
r.h.s. of Eq.|CT] changes according to the DM density field in each 
direction. The set of values {zj f } corresponding to the pixels in a 
given map and their average z re f can be determined by comparing, 
pixel per pixel, a full map of the DM-induced intensity (using fully 
Eq. Q3 and a map containing only the result of the integral in the 
r.h.s. of Eq. lC"Tl ( which we call a /-map). This needs to be done sep- 
arately for all the different shells since z ref changes shell by shell, 
and then combined to produce the total observed emission map. 

There are no approximations made up to this point. We ar- 
gue now that a map for a generic particle physics model can be 

21 An analogous equation for the case of decaying DM can be written, of 
course. 



reconstructed multiplying the J-map by the corresponding particle 
physics factor evaluated at the set of {zL ( ] obtained for our refer- 
ence case as described above. Moreover, to a very good approxi- 
mation, the pixel average value z re f can be used instead of the full 
set {zLf}. We test these arguments by using this technique to recon- 
struct the gamma-ray map for an annihilating DM candidate with 
a mass of 2 TeV, a cross-section of (cr ann v) = 3 x 10~ 26 cm 3 s~' and 
annihilating only into r leptons obtained from a sky map for a ref- 
erence case of a DM candidate with a mass of 200 GeV, the same 
cross section, but with an annihilation channel into bottom quarks. 
We then compare the reconstructed map with one corresponding to 
the same particle physics scenario but obtained directly from the 
complete map-making code. The test is restricted to the EG-MSII 
component for an energy of 10 GeV and to the simulation output 
corresponding to z = 3.1 and Az = 0.25, which is larger than the 
shell with the largest redshift "thickness" we consider in this work. 
We find that the reconstructed map has essentially the same APS 
as the original maps, and the average intensities of the two maps 
agree at the level of 1%. This reconstruction method is not only 
precise when the reconstructed map is obtained accounting for the 
pixel dependence of z re f, but also when the constant average value 
z re f is used for all pixels. 

We are then confident that this procedure can be used to re- 
construct maps of the extragalactic gamma-ray emission for any 
particle physics model. 



APPENDIX D: ANISOTROPY FROM UNRESOLVED 
GALACTIC SUBHALOS 

In the presen t section we d escribed how we implement the method 
described in lAndol d2009h to compute the APS of galactic unre- 
solved subhalos. 

For the subhalo radial distribution we adopt an Einasto profile 
with parameters chosen to match those of the Aq-A-1 main halo: 
M 20 o = 9.4 x 10-" Ms, r_ 2 = 199 kpc, c_ 2 = 1-24, and a = 0.678, 
with the normalization set by the fraction of the smooth halo mass 
M200 in subhalos f m \>. We require / su b = 0.136 for subhalo masses 
in the range 1.7 x 10 s to 10 10 M Q , which is the fraction of the halo 
mass found in resolved subhalos in Aq-A-1; extrapolating the mass 
function to M mm below the minimum resolved subhalo mass leads 
to larger values of / sub . We take the subhalo mass function slope to 
be -1.9, and evaluate the anisotropy for several values of M mm . 

The substructure luminosity function for annihilation is deter- 
mined by assuming the subhalo luminosity is related to the subhalo 
mass by L(M mb ) = AK(M, ub /M e f, with K = b A (av)N y l{2m 2 x ) and 
A a normalization set related to the "astrophysical factor". We con- 
sider two sets of the mass-luminosity parameters (A and /?), chosen 
to reproduce the LOW and HIGH cases in the text. The HIGH case 
extrapolates L(M m b) to M mm using the same relation found to fit the 
resolved subhalos in Aq-A-1; the mass-luminosity relation is cal- 
ibrated to the measured mass-concentration relation and assumes 
each subhalo is well-described by a NFW density profile. For the 
HIGH case we take A = 6.48 x 10 9 M 2 ub kpc~ 3 and p = 0.77. The 
LOW case assumes A = 3.21 x 10 8 M 2 ub kpc" 3 and j3 = 0.86 for 
subhalos with M m b < 1.7 x 10 5 M Q , and the same parameters as the 
HIGH case for subhalos with M sub > 1.7 x 10 5 M Q . The LOW case 
corresponds to a scenario in which subhalo concentrations increase 
more mildly with decreasing subhalo mass, and hence in the LOW 
case the contribution to the intensity and APS from low-mass sub- 
halos is reduced relative to the HIGH case. For decay, the subhalo 
luminosity is always directly proportional to the subhalo mass. 
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We calculated the APS from unresolved subhalos after mask- 
ing the region with \b\ < 30°. We find that for this subhalo model, 
the contribution to the total intensity APS from unresolved subha- 
los for both annihilation and decay is small. For annihilation this 
contribution is ~ 10% of the contribution from the resolved sub- 
halos for the HIGH case, and ~ 5% for the LOW case. For both 
the LOW and HIGH cases the majority of this contribution from 
unresolved subhalos originates from subhalos with masses above 
~ 10 3 M Q . For decay we find that the contribution from unre- 
solved subhalos is at most a few percent of the resolved subhalo 
anisotropy. Since these contributions are small compared to other 
sources of uncertainty in the APS, we do not include them. 
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